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Abstract 

Atomistic simulations are used to test the equations of continuum contact mechanics in nanome- 
ter scale contacts. Nominally spherical tips, made by bending crystals or cutting crystalline or 
amorphous solids, are pressed into a flat, elastic substrate. The normal displacement, contact 
radius, stress distribution, friction and lateral stiffness are examined as a function of load and 
adhesion. The atomic scale roughness present on any tip made of discrete atoms is shown to have 
profound effects on the results. Contact areas, local stresses, and the work of adhesion change 
by factors of two to four, and the friction and lateral stiffness vary by orders of magnitude. The 
microscopic factors responsible for these changes are discussed. The results are also used to test 
methods for analyzing experimental data with continuum theory to determine information, such 
as contact area, that can not be measured directly in nanometer scale contacts. Even when the 
data appear to be fit by continuum theory, extracted quantities can differ substantially from their 
true values. 

PACS numbers: 81.40.Pq 68.35.Np 62.20.Dc 68.37.Ps 
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I. INTRODUCTION 



There has been rapidly growing interest in the behavior of materials at nanometer scales 
One motivation is to construct ever smaller machines P], and a second is to improve 
material properties by controlling their structure at nanometer scales P|. For example, 
decreasing crystallite size may increase yield strength by suppressing dislocation plasticity, 
and material properties may be altered near free interfaces or grain boundaries. To make 
progress, this research area requires experimental tools for characterizing nanoscale prop- 
erties. Theoretical models are also needed both to interpret experiments and to allow new 
ideas to be evaluated. 

One common approach for measuring local properties is to press tips with characteristic 
radii of 10 to 1000 nm into surfaces using an atomic force microscope (AFM) or nanoindenter 
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15|. Mechanical properties are then extracted from the 



measured forces and displacements using classic results from continuum mechanics |16|. A 
potential problem with this approach is that continuum theories make two key assumptions 
that must fail as the size of contacting regions approaches atomic dimensions. One is to 
replace the atomic structure in the bulk of the solid bodies by a continuous medium with 
internal stresses determined by a continuously varying strain field. The second is to model 
interfaces by continuous, differentiable surface heights with interactions depending only on 
the surface separation. Most authors go further and approximate the contacting bodies by 
smooth spheres. 

In a recent paper ^^l, we analyzed the limits of continuum mechanics in describing 
nanometer scale contacts between non-adhesive surfaces with curvature typical of experi- 
mental probes. As in studies of other geometries 



19, 2d, we found that behavior in the 



bulk could be described by continuum mechanics down to lengths as small as two or three 
atomic diameters. However, the atomic structure of surfaces had profound consequences for 
much larger contacts. In particular, atomic-scale changes in the configuration of atoms on 
nominally cylindrical or spherical surfaces produced factor of two changes in the width of 
the contacting region and the stress needed to produce plastic yield, and order of magnitude 
changes in friction and stiffness. 

In this paper we briefly revisit non-adhesive contacts with an emphasis on the role of 
surface roughness. We then extend our atomistic studies to the more common case of 
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adhesive interactions. One important result is that the work of adhesion is very sensitive 
to small changes in the positions of surface atoms. Changes in other quantities generally 
mirror those for non-adhesive tips, and small differences in the magnitude of these effects 
can be understood from geometrical considerations. The results are used to test continuum- 



based methods of analyzing AFM measurements of friction and stiffness 
10, ^|. We show that the models may appear to provide a reasonable description when 



limited information about the true contact structure is available. When the full range 
of information accessible to simulations is examined, one finds that the contact area and 
pressure distributions may be very different than inferred from the models. 

Section ini reviews continuum results for contact without and with adhesion, and briefly 
describes the effect of surface roughness. The methods used in our atomistic simulations 
and the geometries of the tips are described in Sec. IIIII Section |TV| presents results for 
purely repulsive interactions and Sec. |V] describes trends with the strength of adhesion. A 
summary and conclusions are presented in Sec. IVII 

II. CONTINUUM CONTACT MECHANICS 

As noted above, contact mechanics calculations assume that the contacting solids are 
described by continuum elasticity so that the discrete atomic structure can be ignored. In 
most cases the two solids are also assumed to be isotropic with Young's moduli Ei and E2 
and Poisson ratios ui and 1^2 ■ Then the results depend only on an effective modulus E* 
satisfying: 

l/E* = (1 - ul)/E, + (1 - u',)/E2. (1) 

Three-dimensional crystalline solids are not isotropic, but the theories can still be applied 
with an effective E* that depends on orientation and is determined numerically lfi| . 

Continuum theories also neglect the atomic structure of the surface. In most cases the 
surfaces are assumed to be spherical, with radii -Ri and i?2- For elastic, frictionless solids the 
contact of two spheres is equivalent to contact between a sphere of radius R = {R^^ + R2^)~^ 
and a fiat solid ^1^. From Eq. (^, one may then map contact between any two spherical 
surfaces onto contact between a rigid sphere of radius R and a fiat elastic solid of modulus 
E*. This is the case considered in our simulations, and previous results indicate this mapping 
remains approximately correct at atomic scales |l7l |. 



3 



Non-adhesive contact is described by Hertz theory which assumes sohds interact 
with an infinitely sharp and purely repulsive "hard-wall" interaction. The surfaces contact 
in a circular region of radius a that increases with the normal force or load N pushing the 
surfaces together as 

f^V" . (2) 



The normal pressure p within the contact has a simple quadratic dependence on the radial 
distance from the center r: 

, , 2aE* I ^ 

and the surfaces separate slowly outside the contact. The normal displacement of the tip 5 
is related to a by: 

where the subscript H indicates the Hertz prediction and 5h = ^ corresponds to the first 
contact between tip and substrate. 

Adhesion can be treated most simply in the opposite limits of very short-range interac- 
tions considered by Johnson, Kendall and Roberts (JKR) 21 1 and of infinite range interac- 
tions considered by Derjaguin, MuUer and Toporov (DMT) 22*]. The strength of adhesion 
is measured by the work of adhesion per unit area w. In DMT theory the attractive forces 
just produce an extra contribution to the normal force, so that is replaced by + 2tiwR 
in Eqs. ^ and (jl)). JKR theory treats the edge of the contact as a crack tip and calculates 
the stress by adding the crack and Hertz solutions. The normal force in Eq. Q is then 
replaced by A^ + SttwR + [GttwRN + {SirwRyf'^'^ and the equation for 6 is modified (Sec. 



IV B|) . The two approaches lead to very different functional relations between a and A^. For 
example, the contact radius goes to zero at pulloff for DMT theory, but remains finite for 
JKR. They also predict different values of the pulloff force, Nc, where the surfaces separate. 
The normalized pulloff force, Nc/ttwR, is -3/2 in JKR theory and -2 for DMT. Finally, the 
surfaces separate outside the contact with infinite slope in JKR theory, and gradually in 
DMT theory. ^ 

The Maugis-Dugdale (M-D) model j23i] provides a simple interpolation between the JKR 
and DMT limits. The surfaces are assumed to have a hard-wall contact interaction that 
prevents any interpenetration, plus a constant attractive force per unit area, uo, that extends 
over a finite distance ho. The work of adhesion is just the integral of the attractive force. 
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implying ctq/io = w. The M-D model produces coupled equations for the contact pressure 
that can be solved to yield a relation between the load, normal displacement, and area. As 
discussed further in Section |V1 the edge of the contact is broadened by the finite interaction 
range, making it useful to define three characteristic radii that converge to the JKR value 
for a in the limit of short-range interactions. 
Maugis introduced a transition parameter 



A = 



1/3 



(5) 



that measures the ratio of the normal displacement at pulloff from JKR theory to the 



interaction range Hq. Tabor 



2J] had previously defined a similar parameter, fi, that is about 



16% smaller than A for typical interaction potentials [25|. Johnson and Greenwood 25] 
have provided an adhesion map characterizing the range of A over which different models 
are valid. For A > 5 the interaction range is short and JKR theory is accurate, while DMT 
is accurate for A < 0.1. For most materials, both ho and the ratio w/E* are of order 1 nm. 
The JKR limit is only reached by increasing R to macroscopic dimensions of micrometers 
or larger. JKR theory has been tested in experiments with centimeter scale radii using the 
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Scanning probe 



surface force apparatus (SFA) [26] and hemispherical elastomers 
microscope tips typically have R between 10 and 100 nm, and the value of A ~ 0.1 to 1 lies 
between JKR and DMT limits 0] . The same is true in our simulations, where A for adhesive 
tips varies between 0.1 and 0.75. For this reason we will compare our results to M-D theory 
3elow. We also found it useful to use a simple interpolation scheme suggested by Schwarz 
2^. Both he and Carpick et al. jscl have proposed formulae for the contact radius that 
interpolate smoothly between DMT and JKR. These approaches have been attractive in 
analyzing experimental data because of their simple analytic forms. 

No direct measurement of contact area has been possible in nanometer scale single as- 
perity contacts. Instead , th e contact area has been determined by measurements of contact 



stiffness 
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conductance 



0, or friction flQ,^!!^ 
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111. 

' 

The validity of these approaches is not clear [17(, and will be tested below. The stiffness 
against normal displacements of the surfaces can be determined from the derivative of 
with respect to 6 in M-D theory. The tangential stiffness k is normally calculated by assum- 
ing friction prevents sliding at the interface, even though all theories described above assume 
zero friction in calculating the contact area. With this assumption k = 8G*a, where G* is 



the effective bulk shear modulus. Relating the friction to contact area requires assumptions 



about t 



le friction law. Many authors have assumed that the friction is proportional to area 
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|32| and experiments in larger contacts 



but simulations 
33| . I34| | show that this need not be the case. 

The effect of surface roughness on contact has been considered within the continuum 
Q. In general, results ,„„st be obta,ned ...encall. One key parameter is the 
ratio of the root mean squared (rms) roughness of the surface, A, to the normal displacement 
6. When A/ 6 < 0.05, results for nonadhesive contacts lie within a few percent of Hertz 
. As A/6 increases, the contact area broadens and the pressure in the central 
region decreases. Adhesion is more sensitive to roughness js^. The analysis is complicated 
by the fact that A usually depends on the range of lengths over which it is measured. The 
natural upper bound corresponds to the contact diameter and increases with load, while the 
lower bound at atomic scales is unclear. The role of roughness is discussed further in Sees. 
CSIandEl 

III. SIMULATION METHODS 

We consider contact between a rigid spherical tip and an elastic substrate with effective 
modulus E*. As noted above, continuum theory predicts that this problem is equivalent 
to contact between two elastic bodies, and we found this equivalence was fairly accurate 



in previous studies of non-adhesive contact 



17| . To ensure that any deviations from the 



continuum theories described above are associated only with atomic structure, the substrate 
is perfectly elastic. Continuum theories make no assumptions about the nature of the atomic 
structure and interactions within the solids. Thus any geometry and interaction potentials 
can be used to explore the type of deviations from continuum theory that may be produced 
by atomic structure. We use a flat crystalline substrate to minimize surface roughness, and 
use tips with the minimum roughness consistent with atomic structure. The interactions 
are simple pair potentials that are widely used in studies that explore generic behavior js^ . 
They illustrate the type of deviations from continuum theory that may be expected, but 
the magnitude of deviations for real materials will depend on their exact geometry and 
interactions. 

Atoms are placed on sites of a face-centered cubic (fee) crystal with a (001) surface. We 
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define a characteristic length a so that the volume per atom is and the nearest-neighbor 
spacing is 2^/^ a. Nearest-neighbors are coupled by ideal Hookean springs with spring 
constant n. Periodic boundary conditions are applied along the surface of the substrate 
with period L in each direction. The substrate has a finite depth D and the bottom is held 
fixed. For the results shown below, L = 190.5 cr and D = 189.3 a. The continuum theories 
assume a semi-infinite substrate, and we considered smaller and larger L and D to evaluate 
their effect on calculated quantities. Finite-size corrections for the contact radius and lateral 
stiffness are negligible for a/D < 0.1 16j, which covers the relevant range of a/R < 0.2. 
Corrections to the normal displacement are large enough to affect plotted values. We found 



that the leading analytic corrections j37|, |38|, |39| were sufficient to fit our results at large 
loads, as discussed in Sec. IV Bl Note that prev ious simulations of AFM contact have used 
much shallower substrates (D ~ 10 cr) 



hat prev iou 

HHQ, 



This places them in a very 



different limit than continuum theories, although they provide interesting insight into local 
atomic rearrangements. 

Three examples of atomic realizations of spherical tips are shown in Fig. ^ All are 
identical from the continuum perspective, deviating from a perfect sphere by at most a. 
The smoothest one is a slab of f.c.c. crystal bent into a sphere. The amorphous and stepped 
tips were obtained by cutting spheres from a bulk glass or crystal, and are probably more 
typical of AFM tips 2, Q- 

Results for crystaUine tips are very sensitive to the ratio r] 
between their nearest-neighbor spacing and that of the substrate, as well as their crystalline 
alignment 0, 0| • We will contrast results for an aligned commensurate tip with r] = 1 to 
those for an incommensurate tip where rj = 0.94437. To mimic the perfectly smooth surfaces 
assumed in continuum theory, we also show results for a high density tip with t] = 0.05. 
In all cases R = 100 o" ~ 30 nm, which is a typical value for AFM tips. Results for larger 
radius show bigger absolute deviations from continuum predictions, but smaller fractional 
deviations 

Atoms on the tip interact with the top layer of substrate atoms via a truncated Lennard- 
Jones (LJ) potential 0| 

12" 



Vlj = -4ei 



a 

r / \r 



- Vent, r < Tcut (6) 



where characterizes the adhesive binding energy, the potential vanishes for r > rcut, and 
the constant Vent is subtracted so that the potential is continuous at rcut- Purely repulsive 
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interactions are created by truncating the LJ potential at its minimum r^ut = 2^/^ a. Studies 
of adhesion use rcut = 1-5 o or rcut = 2.2 a to explore the effect of the range of the potential. 

In order to compare the effective strength of adhesive interactions and the cohesive energy 
of the solid substrate, we introduce a unit of energy e defined so that the spring constant 
between substrate atoms k = 50 e/cr^. If the solid atoms interacted with a truncated LJ 
potential with e and rcut = 1-5 o", they would have the same equilibrium lattice constant and 
nearly the same spring constant, k = 57 e/cr^, at low temperatures and small deformations. 
Thus ei/e is approximately equal to the ratio of the interfacial binding energy to the cohesive 
binding energy in the substrate. 

The elastic properties of the substrate are not isotropic. We measure an effective modulus 
E* = 55.0 e/a^ for our geometry using Hertzian contact of a high density tip. This is between 
the values calculated from the Young's moduli in different directions. The sound velocity 
is also anisotropic. We find longitudinal sound velocities of 8.5 and 9.5 cr/tij and shear 
velocities of 5.2 and 5.7 cr/tij along the (001) and (111) directions, respectively. Here tij is 
the natural characteristic time unit, t^j = ^Jma'^ je^ where m is the mass of each substrate 
atom. The effective shear modulus for lateral tip displacements is G* = 18.3 e/cx^. 

The simulations were run with the Large-scale Atomic/Molecular Massively Parallel Sim- 



ulator (LAMMPS) code |42ll4S|. The equations of motion were integrated using the velocity- 



Verlet algorithm with time step 0.005 tLj |2^- Temperature, T, only enters continuum 
theory through its effect on constitutive properties, and it is convenient to run simulations 
at low temperatures to minimize fluctuations. A Langevin thermostat was applied to solid 
atoms to maintain T=0.0001 e/ks, where kB is Boltzmann's constant. This is about four 
orders of magnitude below the melting temperature of a Lennard- Jones solid. The damping 
rate was 0.1 t^], and damping was only applied perpendicular to the shding direction in 
friction measurements. 

In simulations, a tip was initially placed about r^ut above the substrate. It was then 
brought into contact by applying a very small load, with the lateral position kept fixed. The 
load was changed in discrete steps and the system was allowed to equilibrate for 350 t^j 
at each load before making measurements. This interval is about 20 times longer than 
the time for sound to propagate across the substrate, and allowed full stress equilibration. 
Results from loading and unloading cycles showed no noticeable hysteresis. To obtain results 
near the pulloff force, we moved the tip away from the substrate at a very slow velocity 
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V = 0.0003 cr/tLj and averaged results over small ranges of displacement. This approach 
was consistent with constant load measurements and allowed us to reach the region that is 
unstable at constant load. 

To compare to continuum predictions we calculated the stresses exerted on the substrate 
by the tip. The force from the tip on each substrate atom was divided by the area per atom 
to get the local stresses. These were then decomposed into a normal stress or pressure p on 
the substrate, and a tangential or shear stress Tgur- The continuum theories described in Sec. 
im assume that the projection of the force normal to the undeformed substrate equals the 
projection normal to the locally deformed substrate. This is valid in the assumed limits of 
a/R << 1 and Tgur = 0. It is also valid for most of our simulations (within < 2%), but not for 
the case of bent commensurate tips where Tgur becomes significant. Normal and tangential 
stresses for bent commensurate tips were obtained using the local surface orientation of the 
nominally spherical tip. Correcting for the orientation changed the normal stress by less 
than 5% of the peak value, and the shear stress by less than 20%. 

Friction forces are known to vary with many parameters [4^. Of particular concern is the 
dependence on extrinsic quantities such as the stiffness of the system that imposes lateral 
motion. Results at constant normal load are often very different than those at fixed height, 
motion perpendicular to the nominal sliding direction can raise or lower friction, and the 
kinetic friction can be almost completely eliminated in very stiff systems 50|, |51|. A full 
re-examination of these effects is beyond the scope of this paper. Our sole goal is to measure 
the friction in a consistent manner that allows us to contrast the load dependent friction for 
different tip geometries and minimizes artifacts from system compliance. 

In friction simulations, the tip is sheared at a constant low velocity v' = 0.01 cr/tLj 
along the (100) direction with a constant normal load. This is typical of AFM experiments 
where the low normal stiffness of the cantilever leads to a nearly constant normal load, and 
the high lateral stiffness limits lateral motion in the direction perpendicular to the sliding 
direction. The measured friction force varies periodically with time as the tip moves by 
a lattice constant of the substrate. The time-averaged or kinetic friction during sliding is 
very sensitive to both lateral stiffness and load jsij. We focus instead on the peak force, 
which is less sensitive. In the limit of low sliding velocities this would correspond to the 
static friction. For bent and stepped commensurate tips there is a single strong friction 
peak. For incommensurate and amorphous tips, there may be multiple peaks of different size 
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corresponding to escape from different metastable energy minima |32l.l5Cl|. Tfie static friction 
was determined from the highest of these friction peaks, since lateral motion would stop at 
any lower force. With a single peak per period, the time between peaks is ~ a /v' = 100 tij- 
This is several times the sound propagation time, and the measured force should be close 
to the static friction. For incommensurate tips the time between peaks was an order of 
magnitude smaller and dynamic effects may be more significant. However, they are not 
expected to affect the load dependence significantly, and are much too small to affect the 
dramatic difference between incommensurate and other tips. 

The total lateral stiffness of the system, k, corresponds to the derivative of F with lateral 
tip displacement evaluated at a potential energy minimum. Since the tip is rigid, k is 
determined by displacements in the substrate and at the interface. The interfacial stiffness 



ki and substrate stiffness fcgub add in series because stress is 
interactions to the substrate. Thus the total stiffness is 17 



transmitted through interfacial 



k-' = Kl + K' ■ (7) 

If the tip were not rigid, it would also contribute a term to the right-hand-side of Eq. ((7|). 

We evaluate k from the derivative of F during sliding, making the assumption that the 
results are in the quasisatic limit. For bent and stepped commensurate tips there is a single 
potential energy minima, and for amorphous tips one minimum dominated the periodic 
force. For incommensurate tips, there are many closely spaced minima and we evaluate k 
from the derivative in the minimum preceding the largest friction peak. Due to the small 
magnitude of forces and short time intervals, the relative errors in these values are as big 
as 50%. To estimate the lateral stiffness in the substrate, fcsub, we fix the relative positions 
of those substrate atoms that lie inside the contact, and move them laterally at a slow 
velocity. The total force between these atoms and the rest of the substrate is measured 
and its derivative with respect to distance gives the lateral stiffness in the substrate. In 
principal, there might also be some relative displacement between atoms in the contact that 
is not captured by this approach, but the results for the substrate stiffness are consistent 
with continuum predictions. 

Values of the adhesion energy per unit area w were obtained for flat, rigid surfaces of 
the same nominal geometry as the tip. For bent crystal tips (Fig. Q), the tip was just 
flattened back into a crystal. For stepped tips, we used an unstepped crystal with the same 
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spacing and interactions. For amorphous tips, an amorphous sohd was cleaved with a flat 
surface rather than a sphere. The resulting surfaces were then brought into contact with 
the substrate and allowed to equilibrate at zero load. At the low temperatures used here, 
the adhesion energy is just the potential energy difference between contacting and separated 
configurations. 

IV. NONADHESIVE CONTACTS 
A. Pressure distribution 

Figure |21 contrasts the distribution of normal pressure p under five tips: (a) dense, (b) 
bent commensurate, (c) bent incommensurate, (d) amorphous and (e) stepped. In each case, 
R = 100 a and the dimensionless load is N/{R^E*) = 0.0018. Hertz theory predicts the 
same pressure distribution (solid lines) for all tips. Points show the actual local pressure on 
each substrate atom as a function of radial distance r from the center of the spherical tip, 
and circles in (c) and (d) show the average over bins of width cr. Clearly, small deviations in 
atomic structure lead to large changes in the mean pressure and the magnitude of pressure 
fluctuations. We find that these deviations become larger as is decreased, and the contact 
radius drops closer to the atomic size. 

One possible source of deviations from Hertz theory is friction, but we find the mean 
tangential forces are small in most cases. The exception is the bent commensurate tip (Fig. 
I2fb)), where the tangential stress rises with r and is comparable to the normal stress near 
the edge of the contact. This result is not surprising given the high friction measured for 
commensurate tips below, and reflects the strong tendency for atoms in the substrate to 
remain locked in epitaxial registry with atoms in the tip. However, the deviation from Hertz 
theory is in the opposite direction from that expected from friction. Since this contact was 
made by gradually increasing the load, friction should decrease the contact size rather than 
broadening it. 

Another possible origin of the deviations from Herts theory is surface roughness. From 
continuum theory (Sec. |H}, this is characterized by the ratio of rms surface roughness A to 
normal displacement 6. The normal displacement for all tips is about the same, 6 ~ 1.5 a, 
but A is difficult to define. The reason is that there is no unique definition of the surface 
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height for a given set of atomic positions. For example, one might conclude that A = for 
the substrate, since all atoms lie on the same plane. However, if a tip atom were moved 
over the surface with a small load, its height would increase as it moved over substrate 
.o,„. and be lowest at sites centered between fout substtate ato^s R. Fot the pa.a™ete.s 
used here, the total height change is about 0.33 a. Similar deviations from a sphere are 
obtained for the bent commensurate and incommensurate tips. The height change decreases 
as the ratio of the nearest-neighbor spacing to the Lennard- Jones diameter for interfacial 
interactions decreases, and is only 0.0007 a for the dense tip. Amorphous and stepped 
tips have additional roughness associated with variations in the position of atomic centers 
relative to a sphere. The total variation is about a, or about three times the height change 
as an atom moves over the substrate. A reasonable estimate is that A/(5 < 0.1 for the bent 
commensurate and incommensurate tips, A/5 < 10~^ for the dense tip, and A/5 ~ 0.3 for 
the amorphous and stepped tips. However, the ambiguity in A is one of the difficulties in 
applying continuum theory in nanoscale contacts. 

The closely spaced atoms on the dense tip approximate a continuous sphere, and the 
resulting pressure distribution is very close to Hertz theory (Fig. Efa)). Results for the 
bent commensurate tip are slightly farther from Hertz theory. The deviations can not be 
attributed to roughness, because fluctuations at a given r are small, and the pressure in the 
central region is not decreased. The main change is to smear the predicted sharp pressure 
drop at the edge of the contact. This can be attributed to the finite range of the repulsive 
potential between surfaces. 

We can estimate the effective interaction range by the change in height of an atom, 
dh = 0.04 a, as p/E* decreases from 0.1 to 0. The effective range is much smaller for the 
dense ti p b ecause ~ 400 times as many atoms contribute to the repulsive potential. In Hertz 
theory [l6|, the separation between surfaces only increases with distance (r — a) from the 
edge of the contact as (8/37r)(r — aY^'^{2aY^'^ / R. Equating this to dh gives r — a 1 a 
for the bent commensurate tip, which is very close to the range over which the edge of the 
contact is shifted from the Hertz prediction. Note that this analysis predicts that the shift 
in the edge of the contact will grow as \/R, and simulations with larger R confirm this. 
However, the fractional change in a decreases as l/\fR. The larger values of pressure at 
low r result from the greater stiffness of the repulsive potential as p increases. All of the 
above effects could be included in continuum theory by changing the form of the repulsive 
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potential 5^ . 

For bent incommensurate and amorphous tips (Fig. 121(c) and (d)), the variations in pres- 
sure at a given r are as large as the mean (circles) [53| . While all atoms on the commensurate 
tip can simultaneously fit between substrate atoms, minimizing fiuctuations in force, atoms 
on the incommensurate tip sample all lateral positions and experience a range of forces at 
a given height. The mean pressure for the incommensurate tip remains fairly close to the 
commensurate results, but there is slightly more smearing at large r due to the variations 
in relative height of substrate and tip atoms. The mean pressure on the amorphous tip 
shows the depression at small r and increase at large r that are expected for rough tips in 
continuum theory The magnitude of the central drop is about 18%, which is consistent 
with A/5 ~ 0.2 in continuum theory (Fig. 13.12 of Ref. [l^). The lack of a noticeable drop 
for incommensurate tips implies that the effective A/6 < .03. The implication is that the 
incoherent height changes on amorphous tips contribute to the effective roughness in con- 
tinuum theory, while the atomic corrugation on bent tips does not. The effective roughness 
in both cases is about 0.1 o" smaller than our estimates for A above. 

Results for stepped tips show the largest deviations from Hertz theory, and they are 
qualitatively different than those produced by random roughness. The terraced geometry of 
this tip (Fig. is closest to that of a flat punch. In continuum theory, the pressure on a 
flat punch is smallest in the center, and diverges as the inverse square root of the distance 
from the edge. The simulation results show qualitatively similar behavior. The main effect 
of atomic structure is to cut off the singularity at a distance corresponding to an atomic 
separation. Similar effects are observed in simulations of other geometries Note 
that the terraces are only flat because the sphere was cut from a crystal that was aligned 



with the substrate. We a. 



so examined tips cut from a crystal that was slightly tilted away 



from the (001) direction 3]. This produces inclined terraces that contact flrst along one 
edge. The resulting pressure distribution is very different, and closest to the continuum 
solution for contact by an asymmetric wedge. 

Figure El has an important general implication about the probability P{p) of flnding a 
local pressure p at a point in the contact. For smoothly curved surfaces, continuum theory 
predicts that the derivative of the pressure diverges at the edge of the contact Thus 
P{p) ^ as p ^ [5^. The flnite resolution at atomic scales always smears out the 
change in p, leading to a non-zero value of -P(O). Indeed, the approximately constant value 
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of dp/dr near the contact edge in Fig. |21 leads to a local maximum in P at j9 = 0. Similar 
behavior is observed for random! 



for piecewise planar surfaces j56|, |57 1 . 



y rough atomic contacts j55|] and in continuum calculations 



Plastic deformation is usually assumed to occur when the deviatoric shear stress exceeds 
the yield stress of the material. In Hertz theory, reaches a maximum value at a depth of 
about 0.5a. The pressure variations at the surface shown in Fig. |211ead to changes in both 



the magnitude and position of the peak shear stress [17[ . Factors of two or more are typical 
for amorphous and stepped tips. Thus tip geometry may have a significant impact on the 
onset of yield. Of course atomistic effects also influence the yield stress at nanometer scales, 
and a full evaluation of this effect is left to future work. Saint- Venant's principal implies 
that the pressure distribution should become independent of tip geometry at depths greater 
than about 3a, but the shear stress at these depths is substantially smaller than peak values 
and yield is unlikely to occur there. 

B. Variations with load 

Figure El shows the load dependence of (a) normal displacement, (b) radius, (c) friction 
and (d) lateral stiffness for the same tips as Fig. |21 Each quantity is raised to a power 
that is chosen so that Hertz theory predicts the indicated straight line. A small finite-depth 
correction (~ 2%) is applied to the Hertz prediction for 6 (Eq. (jH))). 



As also found for cylindrical contacts 



17[, the normal displacement shows the smallest 



deviation from Hertz theory because it represents a mean response of many substrate atoms. 
Results for all bent crystals are nearly indistinguishable from the straight line. Results 
for the stepped surface are lower at small loads. Since the entire tip bottom contacts 
simultaneously, it takes a larger load to push the tip into the substrate. The amorphous 
results are shifted upwards by a fairly constant distance of about 0.2 cr. We define the origin 
of 6 as the tip height where the first substrate atom exerts a repulsive force on the tip. This 
is strongly dependent on the height of the lowest tip atom, while subsequent deformation is 
controlled by the mean tip surface. Agreement with Hertz is greatly improved by shifting 
the amorphous curve by this small height difference. Note that the zero of 6 is difficult to 
determine experimentally and is usually taken as a fit parameter. If this is done, even results 
for the amorphous system yield values of R and E* that are within 10% of their true values. 
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Thus measurements of 5, interpreted with continuum theory for spheres, can provide useful 
values of elastic properties at nanometer scales. 

As expected from the observed pressure distributions (Fig. |2)), the contact radius is 
generally larger than the Hertz prediction. The shift is smallest for the dense tip because it 
approximates a continuous surface and the high density leads to a repulsive potential that 
rises more than a hundred times more rapidly than for other tips. Results for bent crystal 
and amorphous tips are shifted upwards by fairly load-independent offsets of ~ 1 — 3 cr, 
leading to large fractional errors at low loads (up to 100%). The stepped crystal shows 
qualitatively different behavior, with a rising in discrete steps as sequential terraces come 
into contact. Note that the size of the first terrace is not unique, but depends on the registry 
between the bounding sphere and crystalline lattice {3]. Larger deviations may be observed 
when the first step has very few atoms. Such tips may be more likely to be chosen for AFM 
studies because they tend to give sharper images. 

In order to predict the friction between surfaces, one must make assumptions about how 
F depends on area and load. The straight line in Fig. Efc) corresponds to a friction force 
that is proportional to load. Static friction values for bent and stepped commensurate 
surfaces are consistent with this line and a coefficient of friction u = F/N = 0.63. Analytic 

"1 n 

32| and simple numerical |58|| models show that this is a general feature of mated surfaces 
where each tip atom has the same lateral position relative to substrate atoms. The friction 
on amorphous and incommensurate surfaces is always lower and scales more closely with 
the contact indicated by broken line fits to F oc N"^/^ and discussed further in 

Sec. IV Bl [sol] . Many authors have made this assumption in fitting experimental data, but 
it is not obvious why it should hold. The friction per unit area between fiat amorphous 
surfaces decreases as the square root of the area, but rises linearly with the normal pressure 

3^ . Wenning and Miiser have noted that these two factors may combine for spherical 
contacts to give a net friction that rises linearly with area j^^]- However, their argument 
would predict that the frictional force in a cylinder-on-fiat geometry would not scale with 
area, and our previous simulations found that it did [17]. Continuum theory predicts that 
the lateral stiffness k = 8G*a, and should follow the straight line in Fig. Efd). Measured 
values of the total stiffness (open symbols) are always substantially lower. This is because 
continuum theory assumes that there is no displacement at the interface, only within the 
substrate. In reality, the frictional force is always associated with atomic scale displacements 
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of interfacial atoms relative to a local energy minimum j50|, |58[ . The derivative of force with 
displacement corresponds to an interfacial stiffness fcj that adds in series with the substrate 



contribution (Eq. (|7j)) [17|, |5l|. Our numerical results show that k[ can reduce k by more 
than an order of magnitude, particularly for tips where F is small. We also evaluated the 
substrate stiffness ksub by forcing all contacting substrate atoms to move together. These 
results (filled symbols) lie much closer to the Hertz prediction. Only the stepped tip shows 
large discrepancies, and these are correlated with the large deviation between the measured 
and predicted contact radii. 

V. ADHESIVE CONTACTS 
A. Pressure distribution 

Figure ID compares the calculated pressure distribution in adhesive contacts with the 
Maugis-Dugdale prediction (lines). A bent commensurate tip was used to minimize devia- 
tions from continuum predictions for a sphere. Results for two different rcut are presented to 
indicate difficulties in fitting longer-range interactions to M-D theory. The work of adhesion 
was calculated for unbent surfaces (Sec. IIII|) with Cj/e = 0.5, yielding w = 1.05 e/cr^ and 
1.65 e/o"^ for rcut = 1-5 cr and 2.2 a, respectively. This leaves only one fitting parameter in 
M-D theory. For the dashed lines, the width of the attractive interaction ho = w/ao was 
chosen to coincide with the range of the atomic potential. The dotted line shows a fit with 
ho = 0.8 a for rcut = 2.2 a, which gives better values for the pulloff force, but poorer radii 



(Sec. EM- 

In M-D theory, it is common to identify two radii, a and c, with the inner and outer 



23l |25| . For r < a the surfaces are in 



edges of the plateau in the pressure, respectively 
hard-sphere contact, and for a < r < c they are separated and feel the constant attraction. 
The continuously varying interactions between atoms in our simulations lead to qualitatively 
different behavior. There is no sharp transition where the surfaces begin to separate, and 
the attraction shows a smooth rise to a peak, followed by a decay. To facilitate comparison 
to continuum theories, we introduce the three characteristic radii indicated by arrows for 
each rcut- The innermost, r^, corresponds to the point where the interaction changes from 
repulsive to attractive and can be calculated in M-D theory. The outermost, Vc, corresponds 
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to c in M-D theory - the point where interactions end. The middle, Vb, corresponds to the 
peak attraction. It was also studied by Greenwood 5^ for contact of a featureless sphere 
that interacted with the fiat substrate through Lennard- Jones interactions. He found lay 
close to the JKR prediction for contact radius at large loads. All three radii converge in the 
limit of repulsive interactions or the JKR limit of an infinitely narrow interaction range. 

At small radii the atomistic results for p lie above M-D theory, and they drop below near 
r/a = 10. These deviations can be attributed to the increasing stiffness of the repulsion 
between tip and substrate atoms with increasing pressure. Just as in the non-adhesive 
case, the stiffer interactions in the center of the tip lead to bigger pressures for a given 
interpenetration. The change in pressure with separation produces less smearing at the 
edge of the repulsive region of the contact than for the non- adhesive case (Fig. IH^b)), and 
the values of Va are very close to M-D theory. This can be understood from the fact that 
surfaces separate much more rapidly in JKR theory (oc (r — a)^^^) than in Hertz theory 
(oc (r — a)^/^) (Sec. I1V|) . Thus the same height change in the repulsive region corresponds 
to a much smaller change in radius. Of course the finite range of the attractive tail of the 
potential leads to a broad region of adhesive forces out to Vc- The continuous variation 
of p in the attractive tail is only crudely approximated by M-D theory. The difficulty in 
determining the optimum choice for ho increases with the range of interactions, as discussed 
further below. 

Figure IHl shows the effect of tip geometry on pressure distribution. We found that the 
work of adhesion was very sensitive to tip geometry. To compensate for this effect, we varied 
ej (Table 1) to yield the same w for the tips in Fig. El Then all tips should have the same 
pressure distribution in continuum theory. The M-D predictions for p are not shown because 
even the bent commensurate tips produce significantly different results (Fig. HJ. Instead, 
we compare other tips to the bent commensurate tip. 

As for non-adhesive tips, local fluctuations in pressure are small for commensurate tips 
(Fig. Efa) and (d)) and comparable to the mean pressure for incommensurate or amorphous 
tips (Fig. Efb) and (c)). Note however, that the fluctuations become smaller in the adhesive 
regime (large r). This is because the potential varies more slowly with height, so fluctuations 
in the separation between atoms have less effect on p. One consequence is that the outer 
edge of the contact is nearly the same for commensurate and incommensurate tips. The radii 
for the amorphous tip are significantly larger than bent tips, presumably because of a much 
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greater effective roughness. As for nonadhesive tips, the mean pressure on incommensurate 
tips is close to the commensurate results. Adhesion reduces the roughness-induced drop 
in pressure in the central region of the amorphous tip. For the stepped tip, the contact 
radius is dominated by the size of the terraces, but adhesive tails are visible outside the 
edge of each terrace. Only is easily defined for the stepped tips. This increases in a 
nearly stepwise manner, and its load dependence is not shown below. For the amorphous 
and incommensurate tips, radii are determined from the mean pressure at a given radius 
(open circles). Errorbars are less than 0.5 a. 



B. Variations of radius and displacement with load 

Figure Efb) compares the measured contact radii for the tips of Fig. |S1 to M-D theory as 
load is varied. The simulation results for (open symbols) decrease with decreasing load 
as Ta (closed symbols) decreases to zero. All interactions in the contact are then adhesive. 
As Tc continues to drop, the area contributing to adhesion drops, and the load rises back 
toward zero. This regime is not considered in the original M-D theory. The extension to 

= by Johnson and Greenwood (JG) 25| is shown by a dashed line in the figure. Along 



this line the stress in the contact is constant, giving = — ctq Trr^. 

As for non-adhesive tips, the contacts tend to be larger than continuum theory predicts. 
However, the shift for bent tips is smaller than in the non-adhesive case, and the commen- 
surate and incommensurate results are closer, as expected from Fig. [SI Larger deviations 
are observed for the amorphous tip, with radii typically 2 or 3 o" larger than predicted. The 
deviation becomes most pronounced at negative loads, where the amorphous tip remains in 
contact well below the predicted pulloff force. 

Figure El^a) compares the value of rf, to the JKR prediction for contact radius. As found 
by Greenwood (3], the numerical results are close to JKR at large loads, but deviate at 
negative loads because M-D and JKR predict different pulloff forces. Since JKR assumes a 
singular adhesive stress at the radius of the contact, it seems natural that its predictions lie 
closest to the position of the peak tensile stress. 

Figure [3 shows the calculated radii for bent and amorphous tips with the same interaction 
energy ej/e = 0.5. The small changes in tip geometry lead to a roughly four fold variation 
in both w and Nc (Table I). The largest w is obtained for commensurate tips because all 
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atoms can optimize their binding coherently. Atoms on incommensurate tips sample all 
positions relative to the substrate, and can not simultaneously optimize the binding energy. 
The larger height fluctuations on amorphous tips lead to even greater reductions in w. 

In the simulations, the pulloff force, N^.., corresponds to the most negative load where the 
surfaces remain in contact and the rj can be defined. Its normalized value, Nd ttwR, is equal 
to -3/2 in JKR theory, -2 in DMT theory, and lies between these values for M-D theory. 
Table I shows normalized results for various tips. As expected from the good agreement in 
Figs, iniandl?! results for bent tips lie between JKR and DMT values and can be fit with M-D 
theory. The values for stepped and amorphous tips lie outside the bounds of M-D theory. 
This is an important result because pulloff forces are frequently used to measure the work 
of adhesion. Based on continuum theory, one would expect that the uncertainty in such 
measurements is less than the difference between JKR and DMT predictions. Our results 
show factor of two deviations for stepped tips, which may be common in scanning probe 
experiments. Other simulations showed that the stepped tip values are strongly dependent 
on the size of the first terrace, as well as any tilt of the terraces or incommensurability. 

It might seem surprising that the stepped tip has a smaller pulloff force than the bent 
tip, because the entire first terrace can adhere without any elastic energy penalty. However 
this effect is overcome by the greater contact size for bent tips: The radius of the first 
terrace, ~ 6 cr, is smaller than the values of and at pulloff for bent tips. As the 
adhesion is decreased, the predicted contact size at pulloff will drop below rt and the stepped 
tip may then have a larger pulloff force than predicted. This limit can also be reached by 
increasing the width of the first terrace. For a tip that lies entirely within a sphere of radius 
R, r^<R^ — {R — dY ~ 2dR where d is the layer spacing in the crystal. For our geometry 
this corresponds to rt < 12 a, which is about twice the value for our tip. As noted above, 
terraces with even smaller radii may be preferentially selected for imaging studies and will 
lead to lower \Nc\. 

As Tcut increases, it becomes harder for the M-D model to simultaneously fit both the 
radii and the pulloff force. Figure |H1 shows simulation data for a bent commensurate tip 
with Tcut = 2.2 a. Using the value of /iq = 1.2 a (Fig. reproduces all radii fairly 
well at large loads, but gives a substantially different pulloff force, —906 e/a instead of 
—872 e/cr. Decreasing to 0.8 a fits the pulloff force, and improves the fit to at low A^. 
However, the predicted values for at large N are slightly too high and the values for 
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are shifted far below (~ 2 — 3 cr) simulation data. This failure is not surprising given the 
crude approximation for adhesive interactions in M-D theory. As might be expected, the 
best values for the pulloff force are obtained by fitting the region near the peak in the force, 
rather than the weak tail (Fig. E]). 

It should be noted that for bent crystalline and amorphous tips, all of our results for r 
can be fit accurately to M-D theory if E* , w, and R are taken as adjustable parameters. The 
typical magnitude of adjustments is 10 to 30%, which is comparable to typical experimental 
uncertainties for these quantities in nanoscale experiments. Indeed one of the common goals 
of these experiments is to determine any scale dependence in continuum parameters. For 
this reason it would be difficult to test continuum theory in scanning probe experiments 
even if r could be measured directly. 

Experiments can more easily access the variation of normal displacement with load. This 
requires subtracting the height change due to the normal compliance of the machine control- 
ing the tip, which is difficult for standard AFM cantilevers, but possible in stiffer devices 
3y |3. The a..Me .e.o of . . not Exed. bu. be fl.ed to tW. Fi..e 

shows two measures of the normal displacement in our simulations. One, 6tip, corresponds 
to the experimentally accessible tip displacement. We associate the zero of 5tip with the 
point where the first substrate atom exerts a repulsive force. The second, 5sur, is the actual 
depression of the lowest substrate atom on the surface relative to the undeformed surface. 
The two differ because of the interfacial compliance normal to the surface, which is assumed 
to vanish in continuum theory. 

The simulation results for 6 are more sensitive to the finite sample depth D and lateral 
periodicity L than other quantities. To account for this, the predictions of M-D and JKR 



theory are shifted by the leading analytic corrections ina/D |37L l38l l39| 



5 = (5H(l-&-^)+5adhcsion(l-C^-^), (8) 

where 5h is the Hertz prediction (Eq. (jU), h and d are fit parameters, and ^adhesion = 
-^J2wiia/E* for JKR theory and -{2aQ/ E*)^/c^ - a? for M-D theory. We obtained h = 0.8 
from simulations with dense, non-adhesive tips (Fig. 0]) and d = 0.3 from simulations with 
dense, adhesive tips. Results for dense tips are then indistinguishable from the fit lines in 
Figs. ISfa) and M Values of b and d from numerical studies of continuum equations are of 



the same order l37 



mm, 

but to our knowledge these calculations have not considered our 
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geometry where L D. 

With the finite-size corrections, continuum results for 5 he very close to the simulation 
data for bent tips. Agreement is best for dsnr because neither it nor continuum theory 
include the interfacial compliance. The choice of zero for 5tip is not precisely consistent with 
M-D theory. The first repulsive interaction would correspond to the first non-zero r^, which 
occurs at a slightly negative 5 in M-D theory. For the parameters used here, the shift in 5 
is about 0.24 a, which is about twice the discrepancy between 5tip and M-D theory at large 
loads. This implies that the interfacial compliance produces a comparable correction. 

The effects of interfacial compliance are biggest for the most negative values of 5tip. 
Here (5tip decreases monotonically as N increases to zero. In contrast, 5suv is flat and then 
increases back towards zero. In this regime, = and the only interaction comes from the 
attractive tail of the potential. The net adhesive force gradually decreases (A^ increases) 
as the magnitude of the separation between tip and substrate, /S.5 = t^tip — (5sur, increases 
(inset). Most of this change occurs in 5tip, while the displacement of the surface relaxes 
back to zero as the attractive force on it goes to zero. The JG extension to M-D theory 
is expressed in terms of (5tip and provides a good description of its change with N (dashed 
line), even with the assumption of a constant attractive o"o- 

For amorphous tips, both values of 5 are substantially above M-D theory. The shifts 
are bigger than in the non- adhesive case, about 0.35 cr vs. 0.25 cr. This is correlated with 
the larger than predicted puUoff force. Based on the increase in |A^c|, the effective work 
of adhesion appears to be larger by about 30% than that measured for a flat surface. The 
stronger adhesive contribution to the total normal force leads to a larger value of 5. One may 
understand these changes in terms of the effect of surface roughness. Atomic-scale roughness 
on a nominally flat surface prevents many atoms from optimizing their binding energy. As 5 
and the contact area shrink, the long wavelength height fluctuations become irrelevant and 
no longer prevent the few atoms remaining in the contact from adhering. Thus while the 
large load values of 5 can be flt to continuum predictions with the measured w and a simple 
shift in the origin of 5, the small A^ values correspond to a larger work of adhesion. The 
magnitude of the increase (~ 30%) is modest given that the incommensurate tip has about 
twice as large a w as the amorphous tip for the same interaction energy Cj. 

The data for stepped tips are qualitatively different than the others. As for the non- 
adhesive case, 5 is lower than the continuum prediction at large loads, because the flat tip 
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is harder to push into the substrate. The deviation increases rapidly at negative loads, with 
a sharp drop near = where the contact shrinks to the first terrace. As noted above, 
the radius of the first terrace is smaller than the radius at pulloff predicted by continuum 
theory, and the pulloff force is less than half the predicted value. 



C. Friction and lateral stiffness 



Scanning probe microscopes can most easily detect friction forces and the lateral com- 
phance of the tip/substrate system 0,0]. Figure ITUl shows F as a function of load for five 
different adhesive tips. As for non-adhesive tips, tip geometry has a much larger effect on 
F than other quantities, and values for bent commensurate and incommensurate tips differ 
by two orders of magnitude. 

Since the friction was measured at constant load (Sec. 11111) . only values in the stable 
regime, d6tip/dN > could be obtained. Even in this regime, we found that the tip tended 
to detach at A^ > Nc. This was particularly pronounced for commensurate tips. Indeed the 
bent and stepped commensurate tips detached after the first peak in the friction force for 
all the negative loads shown in Fig. [TUT a). At loads closer to the pulloff force, detachment 
occurred at even smaller lateral displacements. 

As noted above, bent commensurate tips have the strongest adhesion energy because all 
atoms can simultaneously optimize their binding. For the same reason, the adhesion energy 
changes rapidly as atoms are displaced laterally away from the optimum position, allowing 
pulloff above the ex pec ted N^.. The extent of the change depends on the sliding direction 
and registry 40|, |49|, |50| . We consider sliding in the (100) direction (Sec. IIII|) . where atoms 
move from points centered between substrate atoms towards points directly over substrate 
atoms. This greatly reduces the binding energy, leading to detachment at less than half the 
pulloff force. Note that changes in binding energy with lateral displacement lead directly 
to a lateral friction force [s^ and the bent and stepped commensurate tips also have the 
highest friction. We suspect that tips with higher friction may generally have a tendency to 
detach farther above during sliding than other tips. 

At the macro scale, friction is usually assumed to be directly proportional to load. All 
the tips have substantial friction at zero load, due to adhesion. The friction force also varies 
non-linearly with A^, showing discrete jumps for the stepped tip, and curvature for the other 
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tips, particularly near Nc- The curvature at small N is reminiscent of the dependence of 
radius on load. Several authors have fit AFM friction data assuming that F scales with 
contact area, and using a continuum model to determine the function of load 



ll|. The dotted lines in Fig. ^Jshow attempts to fit F using the area 



from JKR theory with w adjusted to fit the pulloff force and the proportionality constant 
chosen to fit the friction at high loads. None of the data is well fit by this approach. 
The solid lines show fits to an expression suggested by Schwarz that allows one to 
simply interpolate between DMT and JKR limits as a parameter ti is increased from to 
1. Reasonable fits are possible for non-stepped tips when this extra degree of freedom is 
allowed. Values of the friction per unit within 25% of those determined below from 

the true area (Table II), but the values of w and ti do not correspond to any real microscopic 
parameters. For example, for the three non-stepped tips where a direct measurement gives 
w = 0.46 e/cr^, the fits give w = 0.41 e/o"^, 0.55 e/a^ and 0.94 e/cr^ for commensurate, 
amorphous and incommensurate tips, respectively. The value of ti varies from 0.1 to 0.5 
to 1.5, respectively. Note that the value of ti = 1.5 for incommensurate tips is outside the 
physical range. In this case, and some experiments Q], the friction rises more slowly than 
the JKR prediction, while M-D theory and simpler interpolation schemes [2^ 30 1 always 
give a steeper rise. Such data seem inconsistent with F scaling with area. 

Our simulations allow us to test the relationship between F and area without any fit 
parameters. However, it is not obvious which radius should be used to determine area. 
Figure ITD shows friction plotted against r^, and r^. The stepped tip is not shown, since 
only Tc is easily defined and it increases in one discontinuous jump over the range studied. 
For all other tips the friction is remarkably linear when plotted against any choice of radius 
squared. In contrast, plots of vs. show significant curvature. 

For F to be proportional to area, the curves should also pass through the origin. This 
condition is most closely met by r^, except for the incommensurate case. The idea that 
friction should be proportional to the area where atoms are pushed into repulsion seems 
natural. The other radii include attractive regions where the surfaces may separate far 
enough that the variation of force with lateral displacement is greatly reduced or may even 
change phase. However, in some cases the extrapolated friction remains finite as goes to 
zero and in others it appears to vanish at finite r^. 

Also shown in Fig. ^2 are results for non-adhesive tips (open circles). Values for in- 
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commensurate and amorphous tips were doubled to make the hnearity of the data more 
apparent. Only the bent commensurate data shows significant curvature, primarily at small 
A^. As noted in Sec. IIVI friction is proportional to load for this tip, and the curvature is 
consistent with this and ~ N"^^^. No curvature is visible when adhesion is added, but the 
fact that the linear fit extrapolates to F = at Ta/cr ~ 4 suggests that the linearity might 
break down if data could be obtained at lower loads. 

Results for adhesive and non-adhesive amorphous tips can be fit by lines through the 
origin (not shown), within numerical uncertainty. The same applies for non-adhesive bent 
incommensurate tips, but adding adhesion shifts the intercept to a positive force at = 0. 
It would be interesting to determine whether this extrapolation is valid. Friction can be 
observed with purely adhesive interactions, since the only requirement is that the magnitude 
of the energy varies with lateral displacement However, it may also be that the friction 



on incommensurate tips curves rapidly to zero as ^ 0. Indeed one might expect that for 
very small contacts the tip should behave more like a commensurate tip, leading to a more 
rapid change of F with area or load. The only way to access smaller is to control the tip 
height instead of the normal load. This is known to affect the measured friction force (49l |. 
and most experiments are not stiff enough to access this regime. 

Figure IT^ shows the lateral stiffness as a function of the three characteristic radii. Except 
for the stepped tip (not shown), k rises linearly with each of the radii. As for non-adhesive 
tips, the slope is much smaller than the value 8G* predicted by continuum theory (solid lines) 
because of the interfacial compliance (Eq. (jZj)). The intercept is also generally different from 
the origin, although it comes closest to the origin for r^. As for friction, it seems that the 
repulsive regions produce the dominant contribution to the stiffness. 

Results for non-adhesive tips are also included in Fig. ^1 They also are linear over the 
whole range, and the fits reach A; = at a finite radius a^. This lack of any stiffness between 
contacting surfaces seems surprising. Note however that linear fits to Fig. Efb) also would 
suggest that the radius approached a non-zero value, oq, in the limit of zero load. Moreover, 
the values of oq and follow the same trends with tip geometry and have similar sizes. 

The finite values of Qq and can be understood from the finite range of the repulsive part 
of the interaction. As long as atoms are separated by less than r^ut there is a finite interaction 
and the atoms are considered inside a. However, the force falls rapidly with separation, and 
atoms near this outer limit contribute little to the friction and stiffness. If 6h is the distance 
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the separation must decrease to get a significant interaction, then = {2R5hY/'^ at the 
point where the first significant force is feh. Taking the estimate of 5h = 0.04 a from Sec. 
IIVI then oq ~ 3 cr, which is comparable to the observed values of Oq and a^. The shift is 
smaller for the adhesive case because there are still strong interactions when goes to zero. 
The larger shifts for the amorphous tip may reflect roughness, since the first point to contact 
need not be at the origin. We conclude that the linear fits in Fig. El go to zero at finite 
radius, because overestimates the size of the region that makes significant contributions to 
forces, particularly for non-adhesive tips. Note that the plots for friction with non-adhesive 
tips (Fig. [TT|) are also consistent with an offset, but that the offset appears much smaller 
when plotted as radius squared. 

The slope of the curves in Fig. ^2 can be used to define a differential friction force 
per unit area or yield stress r/ = dF/dTcr^. It is interesting to compare the magnitude of 
these values (Table II) to the bulk yield stress of the substrate Ty. Assuming Lennard- Jones 
interactions, the ideal yield stress of an fee crystal in the same shearing direction is 4 to 10 
ea^, depending on whether the normal load or volume is held fixed. The commensurate tip 
is closest to a continuation of the sample, and the force on all atoms adds coherently. As 
a result tj is of the same order as Ty, even for the non-adhesive tip. Values for adhesive 
amorphous and incommensurate tips are about one and two orders of magnitude smaller 
than Ty, respectively. This reflects the fact that the tip atoms can not optimize their registry 
with the substrate. Removing adhesive interactions reduces Ty by an additional factor of 
about four in both cases. 

In continuum theory, k = 8G*r, and the slope of fits in Fig. could then be used 
to determine the effective shear modulus G*. However, as noted above, the interfacial 
compliance leads to much lower stiffnesses. To illustrate the magnitude of the change we 
quote values of G' = {l/8)dk/dr in Table II. All values are below the true shear modulus 
G* = 18.3 e/cr^ obtained from the substrate compliance alone (Fig. E)). As always, results 
for bent commensurate tips come closest to continuum theory with G'/G* ~ 0.7. Values 
for adhesive amorphous and incommensurate tips are depressed by factors of 3 and 20 
respectively, and removing adhesion suppresses the value for amorphous tips by another 
factor of four. 

Carpick et al. noted that if friction scales with area and k with radius, then the ratio 
F/,' .ho„.d be co„.ta„t gflflQ. Defl„i„g the Metiona. fo.ee per un.t a.ea as and 
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using the expression for k from continuum theory, one finds QAG*"^ F / Tik"^ = Tf. In principal, 
this allows continuum predictions to be checked and Tf to be determined without direct 
measurement of contact size. Figure IT^ shows: 

rf = 64{G*yF/7Te (9) 

as a function of for different tip geometries and interactions. Except for the non-adhesive 
stepped tip, the value of F/k"^ is fairly constant at large loads, within our numerical accuracy. 
Some of the curves rise at small loads because the radius at which F reaches zero in Fig. 
im tends to be smaller than that where k reaches zero in Fig. El These small radii are 
where continuum theory would be expected to be least accurate. Note that the deviations 
are larger for the non-adhesive tips, perhaps because the data extends to smaller radii. 

The data for stepped tips are of particular interest because the contact radius jumps 
in one discrete step from the radius of the first terrace to the radius of the second. The 
friction and stiffness also show discontinuous jumps. Nonetheless, the ratio F/fc^ varies 
rather smoothly and even has numerical values close to those for other tips. The most 
noticeable difference is that the data for the nonadhesive stepped tip rises linearly with 
load, while all other tips tend to a constant at high load. These results clearly demonstrate 
that success at fitting derived quantities like F and k need not imply that the true contact 
area is following continuum theory. 

The curves for r^^ in Fig. El are all much higher than values of the frictional stress Tf 
obtained directly from the friction and area (Table II). Even the trends with tip structure 
are different. The directly measured frictional stress decreases from bent commensurate to 
amorphous to bent incommensurate, while r^^ is largest for the amorphous and smallest 
for the bent commensurate tip. These deviations from the continuum relation are directly 
related to the interfacial compliance ki. The continuum expression for the lateral stiffness 
neglects ki and gives too small a radius at each load. This in turn over-estimates the frictional 
stress by up to two orders of magnitude. Similar effects are likely to occur in experimental 
data. 

Experimental plots of F/ k"^ have been obtained for silicon-nitride tips on mica and sodium 
chloride la LZl SOl and on carbon fibers Isf . Data for carbon fibers and mica in air showed a 

nn 

rapid rise with decreasing at low loads [Si, iSOj • For mica the increase is almost an order of 
magnitude, which is comparable to our results for non-adhesive bent incommensurate tips. 
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This correspondence may seem surprising given that the experiments measured adhesion in 
air. However, the adhesive force was mainly from long-range capillary forces that operate 
outside the area of contact. Following DMT theory, they can be treated as a simple additive 
load that does not affect the contact geometry. In contrast, data for mica in vacuum is 
well fit by JKR theory, implying a strong adhesion within the contact P, [tI . The measured 
value of F/k'^ is nearly constant for this system, just as in our results for most adhesive 
tips. Results for carbon fibers in vacuum ^ show a linear rise like that seen for nonadhesive 
stepped tips. 

From a continuum analysis of the carbon fiber data, the frictiona^stress was estimated 
to be Tf ~ 300 MPa assuming a bulk shear stress of G* = 9.5 GPa [8|. Note that Fig. IT^ 
would suggest Tf/G* ~ 0.1 to 0.3, while the true values (Table II) are as low as 0.0002. The 
data on carbon fibers could be fit with the bulk shear modulus, but data on mica and NaCl 
60^ indicated that G* was 3 to 6 times smaller than bulk values. Our results show that the 
interfacial compliance can easily lead to reductions of this magnitude and a corresponding 
increase in r/, and that care must be taken in interpreting experiments with continuum 
models. 



VI. DISCUSSION AND CONCLUSION 



The results described above show that many different effects can lead to deviations be- 
tween atomistic behavior and continuum theory and quantify how they depend on tip ge- 
ometry for simple interaction potentials (Fig. Q]). In general, the smallest deviations are 
observed for the idealized model of a dense tip whose atoms form a nearly continuous sphere, 
although this tip has nearly zero friction and lateral stiffness. Deviations increase as the 
geometry is varied from a bent commensurate to a bent incommensurate to an amorphous 
tip, and stepped tips exhibit qualitatively different behavior. Tip geometry has the smallest 
effect on the normal displacement and normal stiffness (Fig. El and ^ because they reflect 
an average response of the entire contact. Friction and lateral stiffness are most affected 
(Fig. Elandlini), because they depend on the detailed lateral interlocking of atoms at the 
interface. 

One difference between simulations and continuum theory is that the interface has a finite 
normal compliance. Any realistic interactions lead to a gradual increase in repulsion with 
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separation rather than an ideahzed hard-wall interaction. In our simulations the effective 
range over which interactions increase is only about 4% of the atomic spacing, yet it impacts 
results in several ways. For bent commensurate tips it leads to an increase in pressure in 
the center of the contact (Figs. |21and|^. The pressure at the edge of nonadhesive contacts 
drops linearly over about 2 a, while continuum theory predicts a diverging slope. The width 
of this smearing grows as the square root of the tip radius and leads to qualitative changes in 



the probability distribution of local pressures 
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These effects could be studied 



in continuum theories with soft-wall interactions. The normal interfacial compliance also 
leads to the offset in linear fits of F vs. and k vs. (Figs. ITTl and IT^ . Fits to the friction 
and local stiffness extrapolate to zero at finite values of because atoms at the outer edge 
of the repulsive range contribute to but interact too weakly to contribute substantially 
to F and k. This effect is largest for non-adhesive tips. 

Approximating a spherical surface by discrete atoms necessarily introduces some surface 
roughness. Even bent crystalline tips have atomic scale corrugations, reflecting the variation 
in interaction as tip atoms move from sites nestled between substrate atoms to sites directly 
above. Amorphous and stepped tips have longer wavelength roughness associated with their 
random or layered structures respectively. This longer wavelength roughness has a greater 
effect on the contacts. For non-adhesive interactions, incommensurate and amorphous tips 
have a lower central pressure and wider contact radius than predicted for ideal spheres. 
These changes are qualitatively consistent with continuum calculations for spheres with 
random surface roughness However the effective magnitude of the rms roughness A is 
smaller than expected from the atomic positions. The correlated deviations from a sphere 
on stepped tips, lead to qualitative changes in the pressure distribution on the surface (Fig. 
121 and El- However, these changes are also qualitatively consistent with what continuum 
mechanics would predict for the true tip geometry, which is closer to a fiat punch than a 
sphere. We conclude that the usual approximation of characterizing tips by a single spherical 
radius is likely to lead to substantial errors in calculated properties. Including the true tip 
geometry in continuum calculations would improve their ability to describe nanometer scale 
behavior. Unfortunately this is rarely done, and the atomic-scale tip geometry is rarely 
measured. Recent studies of larger tips and larger scale roughness are an interesting step in 
this direction 

Roughness also has a strong influence on the work of adhesion w (Table I). Values of 



28 



w were determined independently from interactions between nominally flat surfaces. For a 
given interaction strength, commensurate surfaces have the highest w, because each atom can 
optimize its binding simultaneously. The mismatch of lattice constants in incommensurate 
geometries lowers w by a factor of two, and an additional factor of two drop is caused by 
the small (A ~ 0.3 a) height fluctuations on amorphous surfaces. In continuum theory, 
these changes in w should produce nearly proportional changes in puUoff force Nc, and tips 
with the same w and ho should have the same Nc- Measured values of Nc differ from these 
predictions by up to a factor of two. It is particularly signiflcant that the dimensionless 
pulloff force for amorphous and stepped tips lies outside the limits provided by JKR and 
DMT theory. Experimentalists often assume that these bounds place tight limits on errors 
in inferred values of w. 

In the case of amorphous tips the magnitude of Nc is 30% higher than expected. The 
higher than expected adhesion in small contacts may reflect a decrease in effective roughness 
because long-wavelength height fluctuations are suppressed. Stepped tips show even larger 
deviations from continuum theory that are strongly dependent on the size of the flrst terraces 
j^]. Tips selected for imaging are likely to have the smallest terraces and the largest 
deviations from continuum theory. 

Adding adhesion introduces a substantial width to the edge of the contact, ranging from 
the point where interactions flrst become attractive to the outer limits of attractive 
interactions Vc (Fig. E}. As the range of interactions increases, it becomes increasingly 
difficult to fit both these characteristic radii and the pulloff force with the simple M-D 
theory (Fig. IH)). For short-range interactions, good fits are obtained with the measured w 
for bent tips. Data for amorphous tips can only be flt by increasing w, due to the reduction in 
effective roughness mentioned above (Fig. For stepped tips the contact radius increases 
in discrete jumps as successive terraces contact the surface. 

The normal interfacial compliance leads to significant ambiguity in the definition of the 
normal displacement as a function of load (Fig. Continuum theory normally includes 
only the substrate compliance, while experimental measures of the total tip displacement 
5tip include the interfacial compliance. The substrate compliance was isolated by following 
the displacement of substrate atoms, 6sut and found to agree well with theory for bent tips 
in the repulsive regime. Johnson and Greenwood's extension of M-D theory j25| includes the 
interfacial compliance in the attractive tail of the potential. It provides a good description of 
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Stip in the regime where = 0. Here Sup — 4ur increases to the interaction range ho. Results 
for amorphous tips show the greater adhesion noted above. Stepped tips follow continuum 
theory at large loads but are qualitatively different at negative loads. 

The most profound effects of tip geometry are seen in the lateral stiffness k and friction 
F, which vary by one and two orders of magnitude respectively. Continuum theories for k 
do not include the lateral interfacial compliance ki. This adds in series with the substrate 
compliance fcsub (Eq- ©)• Except for commensurate tips, ki « ksnh and the interface 



dominates the total stiffness 



17[ . Experiments have also seen a substantial reduction in the 



expected lateral stiffness from this effect 



The friction on non-adhesive commensurate tips (bent or stepped) increases linearly with 
load, as frequently observed for macroscopic objects. In all other cases, F is a nonlinear 
function of load. Our ability to directly measure contact radii allowed us to show that F 
scales linearly with contact area for incommensurate, amorphous and adhesive bent com- 
mensurate tips. These tips also show a linear scaling of k with radius. While these scalings 
held for any choice of radius, the linear fits are offset from the origin. It appears that the 
effective area contributing to friction and stiffness is often a little smaller than the area of 
repulsive interactions corresponding to r^. As noted above, the offset from appears to 
correspond to the finite range over which repulsive forces rise at the interface. 

ExpeHn,e„ta. data Q ff B fl 3 3 U Q for fric.„„ a„d stiffness have bee„ fit to 
continuum theory with the assumptions that F oc and oc r, but without the offsets 
seen in Figs. ^2 and El We showed that our data for bent and amorphous tips could be fit 
in this way (Fig. fTIUl . but that the fit parameters did not correspond to directly measured 
values. This suggests that care should be taken in interpreting data in this manner. 

We also examined the ratio F/k'^. In continuum theory, this is related to the friction per 
unit area rj^ through Eq. Q. Our results for F/k"^ (Fig. IT^ show the range of behaviors 
observed in experiments, with a relatively constant value for adhesive rapid increase 

at low loads in some nonadhesive cases, and a linear rise for non-adhesive stepped tips. 
The directly measured values of Tf (Table II) are smaller than tJ^ by up to two orders 
of magnitude, and have qualitatively different trends with tip geometry. The difference is 
related to a reduction in the stiffness k due to interfacial compliance. This reduces the 
inferred value of bulk shear modulus G' and increases the calculated contact area at any 
given area. We expect that experimental results for F/ k"^ will produce similar overestimates 
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of the true interfacial shear force. 

It remains unclear why F and k should follow the observed dependence on Tq. Analytic 
arguments for clean, flat surfaces indicate that F is very sensitive to structure, with the forces 
on commensurate, incommensurate and disordered surfaces scaling as different powers of 
I . Only when glassy layers are introduced between the surfaces, does the friction 



area 
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scale in a universal manner 
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64( 1 ■ Wenning and Miiser j3l| have argued that the 



friction on clean, amorphous tips rises linearly with area because of a cancellation of two 
factors, but have not considered k. Naively, one might expect that the length over which 
the force rises is a constant fraction of the lattice spacing and that k is proportional to 
F. However, the friction traces change with load and do not always drop to zero between 
successive peaks. We hope that our results will motivate further analytic studies of this 
problem, and simulations with glassy films and more realistic potentials. 

While we have only considered single asperity contacts in this paper, it is likely that the 
results are relevant more broadly. Many experimental surfaces have random roughness on 
all scales that can be described by self-affine fractal scaling. Continuum models of contact 
between such surfaces show that the radius of most of contacts is comparable to the lower 



length scale cutoff in fractal scaling 5^ |57| . This is typically less than a micrometer, sug- 
gesting that typical contacts have nanometer scale dimensions where the effects considered 
here will be relevant. 
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TABLE I: Relations between interaction strength e^, work of adhesion w, N^. and dimensionless 
puhoff force Nc/{itwR) for tips with different atomic-scale geometries. The last column gives the 
dimensionless pulloff force from M-D theory, N^^~^, with the measured w and /iq = 0.5 a. Values 
are accurate to the last significant digit. 



Tip geometry 


c 




lA^cl(f) 


\Nc\ 


\Ni'-°\ 


Commensurate 


0.213 


0.46 


256 


1.77 


1.74 




0.5 


1.05 


569 


1.72 


1.64 


Incommensurate 


0.535 


0.46 


258 


1.79 


1.74 




0.5 


0.45 


238 


1.68 


1.74 


Amorphous 


1.0 


0.46 


326 


2.26 


1.74 




0.5 


0.23 


136 


1.88 


1.83 


Stepped 


0.213 


0.46 


104 


0.72 


1.74 




0.5 


1.05 


168 


0.51 


1.64 



TABLE IL Frictional stresses tj and apparent moduli G evaluated from the derivatives of fits in 
Figs. 111! and [T^ respectively, as a function of tip geometry and work of adhesion w. Values of G' 
are less than the effective shear modulus of the substrate, G* = 18.3 e/a^, and frictional stresses 
are much smaller than expected from continuum expressions (Fig. I13() . Statistical errorbars are 
indicated in parentheses. 



Tip geometry 


w{ea ^) 


Tfiea 3) 


G'(ec7"3) 


Commensurate 





1.35(7) 


13.7(4) 




0.46 


1.82(4) 


12.8(3) 


Incommensurate 





0.0044(7) 


0.70(5) 




0.46 


0.0151(8) 


1.0(3) 


Amorphous 





0.056(4) 


1.53(5) 




0.46 


0.24(2) 


4.3(6) 
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FIG. 1: Snapshots of atoms near the center regions (diameter 50 a) of spherical tips with average 
radius i?=100 a. Prom top to bottom, tips are made by bending a crystal, cutting a crystal 
or cutting an amorphous solid. Three ratios 77 of the atomic spacing in bent crystals to that 
in the substrate are considered; a dense case r/ = 0.05, a commensurate case r/ = 1, and an 
incommensurate case ry = 0.94437. The step structure of cut crystalline tips is not unique, leading 
to variations in their behavior. 
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FIG. 2: (Color online) Local normal pressure vs. radius for five different tip geometries (a) dense 
tip, (b) bent commensurate crystal, (c) bent incommensurate crystal, (d) amorphous, and (e) 
stepped crystal. All are non-adhesive, and have the same nominal radius R = 100 a, e^/e = 1, and 
normalized load N /{E^E*) = 0.0018. Solid lines show the prediction of Hertz theory, dots show 
the pressure on each surface atom, and circles in (c) and (d) show the mean pressure in radial bins 
of width a. Squares in (b) show the component of the tangential force directed radially from the 
center of the contact. The azimuthal component is nearly zero, and tangential forces are much 
smaller for other tips. 



37 




FIG. 3: Dimensionless plots of powers of (a) normal displacement 5, (b) contact radius a, (c) static 
friction F, and (d) lateral stiffness k vs. the cube root of the normal force N for non-adhesive 
tips with R = 100 a. The powers of 6, a and k are chosen so that continuum theory predicts the 
indicated straight solid lines. In (c), the solid line corresponds to F on N, while dashed and dotted 
lines are fits to F oc iV^/^. The continuum predictions for d and r are followed by a dense tip (stars). 
Also shown are results for bent commensurate (pentagons), bent incommensurate (squares), amor- 
phous (triangles) and stepped commensurate tips (circles). In (d), the total lateral stiffness (open 
symbols) lies well below the continuum prediction because of the interfacial compliance. The stiff- 
ness from the substrate alone (filled symbols) scales with radius, as expected from Hertz theory. 
Numerical uncertainties are comparable to the symbol size. 
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10 20 

r/ G 

FIG. 4: (Color online) Local pressure vs. radius for bent commensurate tips with rcut = 1-5 <t (filled 
squares) and rcut = 2.2 a (open squares) and with R = 100 cr, e^/e = 0.5 and N/{B?E*) = 0.0016. 
Dashed lines show Maugis-Dugdale fits with the measured work of adhesion and /iq fit to the range 
of the interactions. A dotted line shows a fit for rcut = 2.2 a with /iq = 0.8 cr that improves 
agreement with the measured pulloff force (Fig. Characteristic radii r^, r^ and rc are defined 
by the locations where the radially averaged pressure first becomes zero, is most negative, and 
finally vanishes, respectively. 
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FIG. 5: (Color online) Normal pressure vs. radius for four different tip geometries: (a) bent 
commensurate crystal, (b) bent incommensurate crystal, (c) amorphous, and (d) stepped crystal. 
In all cases the tip radius R = 100 a, the normalized load N / {Ft? E*) = 0.0023 and the surface 
energy w = 0.46 e/cj^. Dots show the normal pressure on each surface atom. In (b) and (c), circles 
show the mean pressure in radial bins of width o", and lines show the bent commensurate results 
from (a). Horizontal dashed lines are at zero pressure. 
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n/(r'e*) 

FIG. 6: Variation of contact radii with external load for bent commensurate (squares), bent 
incommensurate (triangles) and amorphous (circles) tips with R = 100 a. In (a), the radius rj, 
where the force is most attractive is compared to JKR theory (solid line) . In (b) , values of ra (filled 
symbols) and Tc (open symbols) arc compared to M-D theory (solid lines) with the JG extension 
for Vc when = (dashed line) . All continuum fits use the independently measured surface energy 
w = 0.46 e/cr^ and an interaction range = 0-5 cr that is consistent with the potential range. 
Numerical uncertainty in radii is comparable to the symbol size. 
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FIG. 7: Contact radii as a function of normal load for (a) bent commensurate, (b) bent incom- 
mensurate, and (c) amorphous tips with = 0.5e. While all have the same interaction strength, 
the work of adhesion varies by more than a factor of four. Broken lines show M-D predictions 
for (triangles) and rc (circles) with the indicated values of w and /iq = 0.5 a. The value of 
(squares) is closest to the JKR result (solid line). Deviations from the continuum predictions show 
the same trends as in Fig. El Numerical uncertainty in radii is comparable to the symbol size. 
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(solid line) to r;, (triangles) for a bent commensurate tip with rcut = 2.2 a. For /iq = 1-2 a (dashed 
line) M-D theory fits both rc and at large loads, but the magnitude of the pulloff force is too 
large. Decreasing /iq to 0.8 a (dotted lines) gives excellent agreement with the pulloff force, but 
not the radii (Fig. 0)). Numerical uncertainty in radii is comparable to the symbol size. 
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FIG. 9: Normal displacement as a function of external load measured from (a) the depression of 
the lowest substrate atom (^sur and (b) the displacement of the tip relative to the height where the 
first substrate atom exerts a repulsive force 5tip- The tips are the same as in Fig. IHl The JKR 
prediction is indicated by dotted lines, the M-D prediction by solid lines, and the JG extension for 
= with dashed lines. All are corrected for the finite dimensions of the substrate as described 
in the text, and uncertainties are comparable to the symbol size. The difference — <^tip — *^sur 
reflects the interfacial compliance, and is shown for the bent incommensurate tip in the inset. 
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FIG. 10: (Color online) Static friction F vs. load N for the indicated tip geometries with same 
systems as in Fig. El Attempts to fit the data by assuming F is proportional to the area predicted 
by JKR theory and the interpolation scheme of Ref. [29| are shown by dotted and solid lines, 
respectively. Dashed lines are separate linear fits for the stepped tip for the cases where one or two 
terraces are in repulsive contact. Numerical uncertainties are comparable to the symbol size. 
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FIG. 11: (Color online) Static friction F plotted against radius squared for the indicated tip 
geometries with same parameters as in Fig. IHI Values of r\ (filled circles), (open squares), and 
(filled triangles) are shown for adhesive tips. Open circles show F vs. c? for non-adhesive tips 
with data in (b) and (c) multiplied by a factor of two for clarity. Dashed lines show unconstrained 
linear fits to each data set. Numerical uncertainties are comparable to the symbol size. 
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FIG. 12: (Color online) Lateral stiffness k plotted against radius for the indicated tip geometries 
with same parameters as in Fig. El Values of (filled circles), r?, (open squares), and Vc (filled 
triangles) are shown for adhesive tips. Open circles show F vs. a for non-adhesive tips. Broken 
lines show unconstrained linear fits to each data set, and the solid line indicates the slope predicted 
by continuum theory. Numerical uncertainties are comparable to the symbol size, except for the 
adhesive data in (b) where they may be up to 50%. 
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FIG. 13: (Color online) Ratio of friction to stiffness squared as a function of load for the indicated 
tip geometries (a) without and (b) with adhesion {w = 0.46 e/o"^). Lines are guides to the eye. 
Numerical uncertainties are comparable to the symbol size, except for the bent incommensurate 
data in (b), where they may be as large as 50%. 
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Atomistic simulations are used to test the equations of continuum contact mechanics in nanometer 
scale contacts. Nominally spherical tips, made by bending crystals or cutting crystalline or amor- 
phous solids, are pressed into a flat, elastic substrate. The normal displacement, contact radius, 
stress distribution, friction and lateral stiffness are examined as a function of load and adhesion. 
The atomic scale roughness present on any tip made of discrete atoms is shown to have profound 
effects on the results. Contact areas, local stresses, and the work of adhesion change by factors of 
two to four, and the friction and lateral stiffness vary by orders of magnitude. The microscopic 
factors responsible for these changes are discussed. The results are also used to test methods for 
analyzing experimental data with continuum theory to determine information, such as contact area, 
that can not be measured directly in nanometer scale contacts. Even when the data appear to be 
fit by continuum theory, extracted quantities can differ substantially from their true values. 
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I. INTRODUCTION 

There has been rapidly growing interest in the behav- 
ior of materials at nanometer scales ij. One motivation 
is to construct ever smaller machines 2] , and a second is 
to improve material properties by controlling their struc- 
ture at nanometer scales 0- For example, decreasing 
crystallite size may increase yield strength by suppress- 
ing dislocation plasticity, and material properties may 
be altered near free interfaces or grain boundaries. To 
make progress, this research area requires experimental 
tools for characterizing nanoscale properties. Theoreti- 
cal models are also needed both to interpret experiments 
and to allovif new ideas to be evaluated. 

One common approach for measuring local properties 
is to press tips with characteristic radii of 10 to 1000 nm 
into surfaces using an atomic force microscope (AFM) or 
nanoindenter |4 ^, «^ 1, 1, iO„ U., JJ, ii, ^ [l| . Me- 
chanical properties are then extracted from the measured 
forces and displacements using classic results from con- 
tinuum mechanics J^]. A potential problem with this 
approach is that continuum theories make two key as- 
sumptions that must fail as the size of contacting re- 
gions approaches atomic dimensions. One is to replace 
the atomic structure in the bulk of the solid bodies by 
a continuous medium with internal stresses determined 
by a continuously varying strain field. The second is 
to model interfaces by continuous, differentiable surface 
heights with interactions depending only on the surface 
separation. Most authors go further and approximate 
the contacting bodies by smooth spheres. 

In a recent paper |17| , we analyzed the limits of con- 
tinuum mechanics in describing nanometer scale contacts 
between non-adhesive surfaces with curvature typical of 
experimental probes. As in studies of other geometries 
we found that behavior in the bulk could 
be described by continuum mechanics down to lengths 



as small as two or three atomic diameters. However, 
the atomic structure of surfaces had profound conse- 
quences for much larger contacts. In particular, atomic- 
scale changes in the configuration of atoms on nominally 
cylindrical or spherical surfaces produced factor of two 
changes in the width of the contacting region and the 
stress needed to produce plastic yield, and order of mag- 
nitude changes in friction and stiffness. 



In this paper we briefly revisit non-adhesive contacts 
with an emphasis on the role of surface roughness. We 
then extend our atomistic studies to the more common 
case of adhesive interactions. One important result is 
that the work of adhesion is very sensitive to small 
changes in the positions of surface atoms. Changes in 
other quantities generally mirror those for non-adhesive 
tips, and small differences in the magnitude of these ef- 
fects can be understood from geometrical considerations. 
The results are used to test continuum-based methods 
of analyzing AFM measurements of friction and stiffness 
mililllElllIillll- We show that the models may 
appear to provide a reasonable description when limited 
information about the true contact structure is available. 
When the full range of information accessible to simu- 
lations is examined, one finds that the contact area and 
pressure distributions may be very different than inferred 
from the models. 



Section Hll reviews continuum results for contact with- 
out and with adhesion, and briefly describes the effect of 
surface roughness. The methods used in our atomistic 
simulations and the geometries of the tips are described 
in Sec. IIIII Section Hvl presents results for purely repul- 
sive interactions and Sec. describes trends with the 
strength of adhesion. A summary and conclusions are 
presented in Sec. IVII 
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II. CONTINUUM CONTACT MECHANICS 

As noted above, contact mechanics calculations assume 
that the contacting solids are described by continuum 
elasticity so that the discrete atomic structure can be 
ignored. In most cases the two solids are also assumed to 
be isotropic with Young's moduli Ei and E2 and Poisson 
ratios vi and 1^2 ■ Then the results depend only on an 
effective modulus E* satisfying: 



l/E* = il-uf)/Ei + il-i^i)/E2. 



(1) 



Three-dimensional crystalline solids are not isotropic, but 
the theories can still be applied with an effective E* that 
depends on orientation and is determined numerically 



Continuum theories also neglect the atomic structure 
of the surface. In most cases the surfaces are assumed 
to be spherical, with radii Ri and R2. For elastic, fric- 
tionless solids the contact of two spheres is equivalent to 
contact between a sphere of radius R = + 
and a flat solid From Eq. QJ, one may then map 
contact between any two spherical surfaces onto contact 
between a rigid sphere of radius R and a flat elastic solid 
of modulus E*. This is the case considered in our simula- 
tions, and previous results indicate this map ping remains 
approximately correct at atomic scales |l7l |. 

Non-adhesive contact is described by Hertz theory , 
which assumes solids interact with an infinitely sharp and 
purely repulsive "hard-wall" interaction. The surfaces 
contact in a circular region of radius a that increases with 
the normal force or load N pushing the surfaces together 
as [3: 



" ^ V 4:E* J 



(2) 



The normal pressure p within the contact has a simple 
quadratic dependence on the radial distance from the 
center r: 



p{r) 



2aE* 
ttR 



1 T 



(3) 



and the surfaces separate slowly outside the contact. The 
normal displacement of the tip 5 is related to a by: 
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a^/R 



(3NR 2/3 
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where the subscript H indicates the Hertz prediction and 
5h ~ ^ corresponds to the first contact between tip and 
substrate. 

Adhesion can be treated most simply in the oppo- 
site limits of very short-range interactions considered by 
Johnson, Kendall and Roberts (JKR) and of infi- 
nite range interactions considered by Derjaguin, MuUer 
and Toporov (DMT) The strength of adhesion is 

measured by the work of adhesion per unit area w. In 



DMT theory the attractive forces just produce an ex- 
tra contribution to the normal force, so that N is re- 
placed by + 2t:wR in Eqs. © and (gj). JKR theory 
treats the edge of the contact as a crack tip and cal- 
culates the stress by adding the crack and Hertz solu- 
tions. The normal force in Eq. ([SJl is then replaced by 

1/2 

N + 3nwR + [6ttwRN + (St twR)'^] and the equation 
for 6 is modified (Sec. IVB|I . The two approaches lead 
to very different functional relations between a and N. 
For example, the contact radius goes to zero at puUoff 
for DMT theory, but remains finite for JKR. They also 
predict different values of the puUoff force, N^, where the 
surfaces separate. The normalized puUoff force, Nc/ttwR, 
is -3/2 in JKR theory and -2 for DMT. Finally, the sur- 
faces separate outside the contact with infinite slope in 
JKR theory, and gradually in DMT theory. 

The Maugis-Dugdale (M-D) model 23] provides a sim- 
ple interpolation between the JKR and DMT limits. The 
surfaces are assumed to have a hard-wall contact interac- 
tion that prevents any interpenetration, plus a constant 
attractive force per unit area, ao, that extends over a 
finite distance Hq. The work of adhesion is just the in- 
tegral of the attractive force, implying cto/iq = w. The 
M-D model produces coupled equations for the contact 
pressure that can be solved to yield a relation between 
the load, normal displacement, and area. As discussed 
further in Section^ the edge of the contact is broadened 
by the finite interaction range, making it useful to define 
three characteristic radii that converge to the JKR value 
for a in the limit of short-range interactions. 

Maugis introduced a transition parameter [23| 



9Rw^ 

2TTE*'^hl 



1/3 



(5) 



that measures the ratio of the normal displacement at 
puUoff from JKR theory to the interaction range /iq. Ta- 
bor 24] had previously defined a similar parameter, /i, 
that is about 16% smaller than A for typical interaction 
potentials Johnson and Greenwood 25] have pro- 

vided an adhesion map characterizing the range of A over 
which different models are valid. For A > 5 the inter- 
action range is short and JKR theory is accurate, while 
DMT is accurate for A < 0.1. For most materials, both 
and the ratio w/E* are of order 1 nm. The JKR limit is 
only reached by increasing R to macroscopic dimensions 
of micrometers or larger. JKR theory has been tested in 
experiments with centimeter scale radii using the surface 
force apparatus (SEA) "2^ and hemispherical elastomers 
21s 2^- Scanning probe microscope tips typically have R 
between 10 and 100 nm, and the value of A 0.1 to 1 lies 
between JKR and DMT limits 0. The same is true in 
our simulations, where A for adhesive tips varies between 
0.1 and 0.75. For this reason we will compare our results 
to M-D theory below. We also found it useful to use a 
simple interpolation scheme suggested by Schwarz j29j . 
Both he and Carpick et al. jSOj have proposed formulae 
for the contact radius that interpolate smoothly between 
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DMT and JKR. These approaches have been attractive 
in analyzing experimental data because of their simple 
analytic forms. 

No direct measurement of contact area has been possi- 
ble in nanometer scale single asperity contacts. Instead, 
the contact area has been determined b y rn easurements 
of contact stiffness i i i |1 11 i EMJllJll IH, 
conductance H, or friction BMM S IpLM El The 
validity of these approaches is not clear [13 , and will be 
tested below. The stiffness against normal displacements 
of the surfaces can be determined from the derivative 
of N with respect to S in M-D theory. The tangential 
stiffness k is normally calculated by assuming friction 
prevents sliding at the interface, even though all theories 
described above assume zero friction in calculating the 
contact area. With this assumption k = 8G*a, where G* 
is the effective bulk shear modulus. Relating the friction 
to contact area requires assumptions about the friction 
law. Many authors have assumed that the friction is pro- 
portional to area 0, IE IS S IE S ^> ^| j but simulations 
[ItL I31I Is^ and experiments in larger contacts IIE 
show that this need not be the case. 

The effect of surface roughness on contact has been 
considered within the continuum framework |l6l | . In gen- 
eral, results must be obtained numerically. One key 
parameter is the ratio of the root mean squared (rms) 
roughness of the surface. A, to the normal displacement 
S. When A/S < 0.05, results for nonadhesive contacts 
lie within a few percent of Hertz theory As A/5 

increases, the contact area broadens and the pressure in 
the central re gion decreases. Adhesion is more sensitive 
to roughness '3^. The analysis is complicated by the 
fact that A usually depends on the range of lengths over 
which it is measured. The natural upper bound corre- 
sponds to the contact diameter and increases with load, 
while the lower bound at atomic scales is unclear. The 
role of roughness is discussed further in Sees. llVl andlVl 



III. SIMULATION METHODS 

We consider contact between a rigid spherical tip and 
an elastic substrate with effective modulus E* . As noted 
above, continuum theory predicts that this problem is 
equivalent to contact between two elastic bodies, and 
we found this equivalence was fa irly accurate in previous 
studies of non-adhesive contact [i3 . To ensure that any 
deviations from the continuum theories described above 
are associated only with atomic structure, the substrate 
is perfectly elastic. Continuum theories make no assump- 
tions about the nature of the atomic structure and inter- 
actions within the solids. Thus any geometry and in- 
teraction potentials can be used to explore the type of 
deviations from continuum theory that may be produced 
by atomic structure. We use a flat crystalline substrate 
to minimize surface roughness, and use tips with the min- 
imum roughness consistent with atomic structure. The 
interactions are simple pair potentials that are widely 



used in studies that explore generic behavior [Sq . They 
illustrate the type of deviations from continuum theory 
that may be expected, but the magnitude of deviations 
for real materials will depend on their exact geometry 
and interactions. 

Atoms are placed on sites of a face-centered cubic (fee) 
crystal with a (001) surface. We define a characteris- 
tic length (T so that the volume per atom is and the 
nearest-neighbor spacing is 2^/^ a. Nearest-neighbors are 
coupled by ideal Hookean springs with spring constant 
K. Periodic boundary conditions are applied along the 
surface of the substrate with period L in each direction. 
The substrate has a finite depth D and the bottom is 
held fixed. For the results shown below, L = 190.5 a and 
D = 189.3 a. The continuum theories assume a semi- 
infinite substrate, and we considered smaller and larger 
L and D to evaluate their effect on calculated quanti- 
ties. Finite-size corrections for the contact radius and 
lateral stiffness are neglig ible for a/D < 0.1 El, which 
covers the relevant range of a/R < 0.2. Corrections to 
the normal displacement are large enough to affect plot- 
ted values. We found that the leading analytic correc- 
tions IIE HE 113 were sufficient to fit our results at large 
loads, as discussed in Sec. IV Bl Note that previous sim- 
ulations of AFM contact have used much shallower sub- 
strates (D ~ 10 cr) HEIHIillillillil. This places 
them in a very different limit than continuum theories, al- 
though they provide interesting insight into local atomic 
rearrangements. 

Three examples of atomic realizations of spherical tips 
are shown in Fig. ^ All are identical from the continuum 
perspective, deviating from a perfect sphere by at most 
a. The smoothest one is a slab of f.c.c. crystal bent 
into a sphere. The amorphous and stepped tips were 
obtained by cutting spheres from a bulk glass or crys- 
tal, and are probably more typical of AFM tips 0, Eg ■ 
Results for crystalline tips are very sensitive to the ra- 
tio r/ between their nearest-neighbor spacing and that 
of the substrate, as well as their crystalline alignment 
[31I Is^ l . We will contrast results for an aligned commen- 
surate tip with 77 = 1 to those for an incommensurate 
tip where rj = 0.94437. To mimic the perfectly smooth 
surfaces assumed in continuum theory, we also show re- 
sults for a high density tip with 77 = 0.05. In all cases 
R — 100 (7 ~ 30 nm, which is a typical value for AFM 
tips. Results for larger radius show bigger absolute devi- 
ations from continuum predictions, but smaller fractional 
deviations 

Atoms on the tip interact with the top layer of sub- 
strate atoms via a truncated Lennard- Jones (LJ) poten- 
tial m 



Vlj = -4e, 



(7)' -(7)' 



-Vr, 



r < Tc 



(6) 



where characterizes the adhesive binding energy, the 
potential vanishes for r > Tcut, and the constant Vcut 
is subtracted so that the potential is continuous at r^ut- 
Purely repulsive interactions are created by truncating 
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FIG. 1: Snapshots of atoms near the center regions (diameter 
50 a) of spherical tips with average radius i?=100 a. From 
top to bottom, tips are made by bending a crystal, cutting a 
crystal or cutting an amorphous solid. Three ratios rj of the 
atomic spacing in bent crystals to that in the substrate are 
considered; a dense case rj = 0.05, a commensurate case rj — 1, 
and an incommensurate case rj = 0.94437. The step structure 
of cut crystalline tips is not unique, leading to variations in 
their behavior. 



the LJ potential at its minimum rcut = 2^/^ a. Studies 
of adhesion use rcut = 1-5 a or rcut = 2.2 a to explore 
the effect of the range of the potential. 

In order to compare the effective strength of adhesive 
interactions and the cohesive energy of the solid sub- 
strate, we introduce a unit of energy e defined so that the 
spring constant between substrate atoms k = 50 e/cr^. If 
the solid atoms interacted with a truncated LJ poten- 
tial with e and rcut — 1-5 cr, they would have the same 
equilibrium lattice constant and nearly the same spring 
constant, k — 57 e/cr^, at low temperatures and small de- 
formations. Thus Ei/e is approximately equal to the ratio 
of the interfacial binding energy to the cohesive binding 
energy in the substrate. 

The elastic properties of the substrate are not isotropic. 
We measure an effective modulus E* = 55.0 e/a^ for 
our geometry using Hertzian contact of a high density 
tip. This is between the values calculated from the 
Young's moduli in different directions. The sound ve- 
locity is also anisotropic. We find longitudinal sound 
velocities of 8.5 and 9.5 cr/tLj and shear velocities of 5.2 
and 5.7 cr/tLj along the (001) and (111) directions, re- 
spectively. Here t^j is the natural characteristic time 
unit, t^j — y/ma^Je, where m is the mass of each sub- 
strate atom. The effective shear modulus for lateral tip 
displacements is G* — 18.3 e/a^. 

The simulations were run with the Large-scale 
Atomic/Molecular Massively Parallel Simulator 
(LAMMPS) code [H El. The equations of mo- 
tion were integrated usin g th e velocity- Verlet algorithm 
with time step 0.005 t^j Temperature, T, only en- 
ters continuum theory through its effect on constitutive 
properties, and it is convenient to run simulations at 



low temperatures to minimize fluctuations. A Langevin 
thermostat was applied to solid atoms to maintain 
T=0.0001 e/ks, where fc^ is Boltzmann's constant. 
This is about four orders of magnitude below the 
melting temperature of a Lennard- Jones solid. The 
damping rate was 0.1 tj^\, and damping was only 
applied perpendicular to the sliding direction in friction 
measurements. 

In simulations, a tip was initially placed about rcut 
above the substrate. It was then brought into contact by 
applying a very small load, with the lateral position kept 
fixed. The load was changed in discrete steps and the sys- 
tem was allowed to equilibrate for 350 t^j at each load 
before making measurements. This interval is about 20 
times longer than the time for sound to propagate across 
the substrate, and allowed full stress equilibration. Re- 
sults from loading and unloading cycles showed no no- 
ticeable hysteresis. To obtain results near the puUoff 
force, we moved the tip away from the substrate at a 
very slow velocity v = 0.0003 cr/tLj and averaged results 
over small ranges of displacement. This approach was 
consistent with constant load measurements and allowed 
us to reach the region that is unstable at constant load. 

To compare to continuum predictions we calculated the 
stresses exerted on the substrate by the tip. The force 
from the tip on each substrate atom was divided by the 
area per atom to get the local stresses. These were then 
decomposed into a normal stress or pressure p on the 
substrate, and a tangential or shear stress Tgui-. The con- 
tinuum theories described in Sec. ITU assume that the pro- 
jection of the force normal to the undeformed substrate 
equals the projection normal to the locally deformed sub- 
strate. This is valid in the assumed limits of a/R « 1 
and Tsur = 0. It is also valid for most of our simulations 
(within < 2%), but not for the case of bent commen- 
surate tips where Tgur becomes significant. Normal and 
tangential stresses for bent commensurate tips were ob- 
tained using the local surface orientation of the nominally 
spherical tip. Correcting for the orientation changed the 
normal stress by less than 5% of the peak value, and the 
shear stress by less than 20%. 

Friction forces are known to vary with many parame- 
ters [49I I . Of particular concern is the dependence on ex- 
trinsic quantities such as the stiffness of the system that 
imposes lateral motion. Results at constant normal load 
are often very different than those at fixed height, motion 
perpendicular to the nominal sliding direction can raise 
or lower friction, and the kinetic friction can be almost 
completely eliminated in very stiff systems 15(1 l5ll |. A 
full re-examination of these effects is beyond the scope 
of this paper. Our sole goal is to measure the friction in 
a consistent manner that allows us to contrast the load 
dependent friction for different tip geometries and mini- 
mizes artifacts from system compliance. 

In friction simulations, the tip is sheared at a constant 
low velocity v' = 0.01 (jjtLj along the (100) direction 
with a constant normal load. This is typical of AFM 
experiments where the low normal stiffness of the can- 
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tilever leads to a nearly constant normal load, and the 
high lateral stiffness limits lateral motion in the direction 
perpendicular to the sliding direction. The measured fric- 
tion force varies periodically with time as the tip moves 
by a lattice constant of the substrate. The time-averaged 
or kinetic friction during sliding is very sensitive to both 
lateral stiffness and load |51|. We focus instead on the 
peak force, which is less sensitive. In the limit of low slid- 
ing velocities this would correspond to the static friction. 
For bent and stepped commensurate tips there is a sin- 
gle strong friction peak. For incommensurate and amor- 
phous tips, there may be multiple peaks of different size 
corresponding to escape from different metastable energy 
minima [sS, Ejl • The static friction was determined from 
the highest of these friction peaks, since lateral motion 
would stop at any lower force. With a single peak per 
period, the time between peaks is ~ a/v' — 100 t^j. 
This is several times the sound propagation time, and 
the measured force should be close to the static friction. 
For incommensurate tips the time between peaks was an 
order of magnitude smaller and dynamic effects may be 
more significant. However, they are not expected to af- 
fect the load dependence significantly, and are much too 
small to affect the dramatic difference between incom- 
mensurate and other tips. 

The total lateral stiffness of the system, fc, corresponds 
to the derivative of F with lateral tip displacement eval- 
uated at a potential energy minimum. Since the tip is 
rigid, k is determined by displacements in the substrate 
and at the interface. The interfacial stiffness fei and sub- 
strate stiffness fcsub add in series because stress is trans- 
mitted through interfacial interactions to the substrate. 
Thus the total stiffness is [itIIsH : 

k-'-Kl + K' ■ (7) 

If the tip were not rigid, it would also contribute a term 
to the right-hand-side of Eq. O . 

We evaluate k from the derivative of F during sliding, 
making the assumption that the results are in the qua- 
sisatic limit. For bent and stepped commensurate tips 
there is a single potential energy minima, and for amor- 
phous tips one minimum dominated the periodic force. 
For incommensurate tips, there are many closely spaced 
minima and we evaluate k from the derivative in the min- 
imum preceding the largest friction peak. Due to the 
small magnitude of forces and short time intervals, the 
relative errors in these values are as big as 50%. To esti- 
mate the lateral stiffness in the substrate, fcsub, we fix the 
relative positions of those substrate atoms that lie inside 
the contact, and move them laterally at a slow velocity. 
The total force between these atoms and the rest of the 
substrate is measured and its derivative with respect to 
distance gives the lateral stiffness in the substrate. In 
principal, there might also be some relative displacement 
between atoms in the contact that is not captured by this 
approach, but the results for the substrate stiffness are 
consistent with continuum predictions. 

Values of the adhesion energy per unit area w were 



obtained for flat, rigid surfaces of the same nominal ge- 
ometry as the tip. For bent crystal tips (Fig. QJ, the tip 
was just flattened back into a crystal. For stepped tips, 
we used an unstepped crystal with the same spacing and 
interactions. For amorphous tips, an amorphous solid 
was cleaved with a flat surface rather than a sphere. The 
resulting surfaces were then brought into contact with 
the substrate and allowed to equilibrate at zero load. At 
the low temperatures used here, the adhesion energy is 
just the potential energy difference between contacting 
and separated configurations. 



IV. NONADHESIVE CONTACTS 
A. Pressure distribution 

Figure 121 contrasts the distribution of normal pressure 
p under five tips: (a) dense, (b) bent commensurate, (c) 
bent incommensurate, (d) amorphous and (c) stepped. 
In each case, i? = 100 cr and the dimensionless load is 
N/{R'^E*) = 0.0018. Hertz theory predicts the same 
pressure distribution (solid lines) for all tips. Points 
show the actual local pressure on each substrate atom 
as a function of radial distance r from the center of the 
spherical tip, and circles in (c) and (d) show the aver- 
age over bins of width a. Clearly, small deviations in 
atomic structure lead to large changes in the mean pres- 
sure and the magnitude of pressure fluctuations. We flnd 
that these deviations become larger as N is decreased, 
and the contact radius drops closer to the atomic size. 

One possible source of deviations from Hertz theory 
is friction, but we find the mean tangential forces are 
small in most cases. The exception is the bent commen- 
surate tip (Fig. |2tb)), where the tangential stress rises 
with r and is comparable to the normal stress near the 
edge of the contact. This result is not surprising given 
the high friction measured for commensurate tips below, 
and reflects the strong tendency for atoms in the sub- 
strate to remain locked in epitaxial registry with atoms 
in the tip. However, the deviation from Hertz theory is in 
the opposite direction from that expected from friction. 
Since this contact was made by gradually increasing the 
load, friction should decrease the contact size rather than 
broadening it. 

Another possible origin of the deviations from Herts 
theory is surface roughness. From continuum theory 
(Sec. this is characterized by the ratio of rms sur- 
face roughness A to normal displacement 5. The normal 
displacement for all tips is about the same, J « 1.5 cr, 
but A is difflcult to define. The reason is that there is no 
unique definition of the surface height for a given set of 
atomic positions. For example, one might conclude that 
A = for the substrate, since all atoms lie on the same 
plane. However, if a tip atom were moved over the surface 
with a small load, its height would increase as it moved 
over substrate atoms and be lowest at sites centered be- 
tween four substrate atoms [50j | . For the parameters used 



6 




5 10 15 

v/a 

FIG. 2: (Color online) Local normal pressure vs. radius for 
five diflterent tip geometries (a) dense tip, (b) bent commensu- 
rate crystal, (c) bent incommensurate crystal, (d) amorphous, 
and (e) stepped crystal. All are non- adhesive, and have the 
same nominal radius -R = 100 (t, e^/e = 1, and normalized 
load NI{R^E*) = 0.0018. Solid lines show the prediction of 
Hertz theory, dots show the pressure on each surface atom, 
and circles in (c) and (d) show the mean pressure in radial 
bins of width a. Squares in (b) show the component of the 
tangential force directed radially from the center of the con- 
tact. The azimuthal component is nearly zero, and tangential 
forces are much smaller for other tips. 



here, the total height change is about 0.33 cr. Similar 
deviations from a sphere are obtained for the bent com- 
mensurate and incommensurate tips. The height change 
decreases as the ratio of the nearest-neighbor spacing to 
the Lennard- Jones diameter for interfacial interactions 
decreases, and is only 0.0007 a for the dense tip. Amor- 
phous and stepped tips have additional roughness as- 
sociated with variations in the position of atomic cen- 
ters relative to a sphere. The total variation is about 
cr, or about three times the height change as an atom 
moves over the substrate. A reasonable estimate is that 
A/5 < 0.1 for the bent commensurate and incommensu- 
rate tips, A/(5 < 10^"^ for the dense tip, and A/(5 ^ 0.3 
for the amorphous and stepped tips. However, the ambi- 
guity in A is one of the difhculties in applying continuum 
theory in nanoscale contacts. 

The closely spaced atoms on the dense tip approximate 
a continuous sphere, and the resulting pressure distribu- 



tion is very close to Hertz theory (Fig. Efa)). Results for 
the bent commensurate tip are slightly farther from Hertz 
theory. The deviations can not be attributed to rough- 
ness, because fluctuations at a given r are small, and the 
pressure in the central region is not decreased. The main 
change is to smear the predicted sharp pressure drop at 
the edge of the contact. This can be attributed to the 
finite range of the repulsive potential between surfaces. 

We can estimate the effective interaction range by the 
change in height of an atom, dh = 0.04 cr, as p/E* 
decreases from 0.1 to 0. The effective range is much 
smaller for the dense tip because ~ 400 times as many 
atoms contribute to the repulsive potential. In Hertz the- 
ory |l6l |. the separation between surfaces only increases 
with distance (r — a) from the edge of the contact as 
(8/37r)(r - afl'^{2a)'^l'^/R. Equating this to dh gives 
r — a \ a for the bent commensurate tip, which is very 
close to the range over which the edge of the contact is 
shifted from the Hertz prediction. Note that this analysis 
predicts that the shift in the edge of the contact will grow 
as and simulations with larger R confirm this. How- 
ever, the fractional change in a decreases as l/Vi?. The 
larger values of pressure at low r result from the greater 
stiffness of the repulsive potential as p increases. All of 
the above effects could be included in continuum theory 
by changing the form of the repulsive potential [52| . 

For bent incommensurate and amorphous tips (Fig. 
12 (c) and (d)), the variations in pres sure at a given r 
are as large as the mean (circles) [S^. While all atoms 
on the commensurate tip can simultaneously fit between 
substrate atoms, minimizing fluctuations in force, atoms 
on the incommensurate tip sample all lateral positions 
and experience a range of forces at a given height. The 
mean pressure for the incommensurate tip remains fairly 
close to the commensurate results, but there is slightly 
more smearing at large r due to the variations in relative 
height of substrate and tip atoms. The mean pressure 
on the amorphous tip shows the depression at small r 
and increase at large r that are expected for rough tips 
in continuum theory |l6j . The magnitude of the central 
drop is about 18%, which is consistent with A/5 ^ 0.2 in 
continuum theory (Fig. 13.12 of Ref. [H). The lack of 
a noticeable drop for incommensurate tips implies that 
the effective A/J < .03. The implication is that the 
incoherent height changes on amorphous tips contribute 
to the effective roughness in continuum theory, while the 
atomic corrugation on bent tips does not. The effective 
roughness in both cases is about 0.1 ct smaller than our 
estimates for A above. 

Results for stepped tips show the largest deviations 
from Hertz theory, and they are qualitatively different 
than those produced by random roughness. The ter- 
raced geometry of this tip (Fig. is closest to that 
of a flat punch. In continuum theory, the pressure on 
a flat punch is smallest in the center, and diverges as 
the inverse square root of the distance from the edge. 
The simulation results show qualitatively similar behav- 
ior. The main effect of atomic structure is to cut off 
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the singularity at a distance corresponding to an atomic 
separation. Similar effects are observed in simulations of 
other geometries 0, 0, |23| ■ Note that the terraces are 
only flat because the sphere was cut from a crystal that 
was aligned with the substrate. We also examined tips 
cut from a crystal that was slightly tilted away from the 
(001) direction '4^. This produces inclined terraces that 
contact first along one edge. The resulting pressure dis- 
tribution is very different, and closest to the continuum 
solution for contact by an asymmetric wedge. 

Figure [21 has an important general implication about 
the probability Pip) of finding a local pressure p at a 
point in the contact. For smoothly curved surfaces, con- 
tinuum theory predicts that the derivative of the pressure 
diverges at the edge of the contact [l^. Thus P{p) 
as p ^ The finite resolution at atomic scales al- 
ways smears out the change in p, leading to a non-zero 
value of P(0). Indeed, the approximately constant value 
of dp/dr near the contact edge in Fig. |21 leads to a local 
maximum in P at p = 0. Similar behavior is observed for 
randomly rough atomic contacts IsBj and in continuum 
calculations for piecewise planar surfaces [s^ [sj ■ 

Plastic deformation is usually assumed to occur when 
the deviatoric shear stress Ts exceeds the yield stress of 
the material. In Hertz theory, Ts reaches a maximum 
value at a depth of about 0.5a. The pressure variations 
at the surface shown in Fig. El lead to changes in both 
the magnitude and position of the peak shear stress ■ 
Factors of two or more are typical for amorphous and 
stepped tips. Thus tip geometry may have a significant 
impact on the onset of yield. Of course atomistic ef- 
fects also influence the yield stress at nanometer scales, 
and a full evaluation of this effect is left to future work. 
Saint- Venant's principal implies that the pressure distri- 
bution should become independent of tip geometry at 
depths greater than about 3a, but the shear stress at 
these depths is substantially smaller than peak values 
and yield is unlikely to occur there. 
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FIG. 3: Dimensionless plots of powers of (a) normal dis- 
placement 5, (b) contact radius a, (c) static friction F, and 
(d) lateral stiffness k vs. the cube root of the normal force 
N for non-adhesive tips with R = 100 a. The powers of 5, a 
and k are chosen so that continuum theory predicts the indi- 
cated straight solid lines. In (c), the solid line corresponds to 
F oc N, while dashed and dotted lines are fits to _F oc N^^^ . 
The continuum predictions for S and r are followed by a dense 
tip (stars). Also shown are results for bent commensurate 
(pentagons), bent incommensurate (squares), amorphous (tri- 
angles) and stepped commensurate tips (circles). In (d), the 
total lateral stiffness (open symbols) lies well below the con- 
tinuum prediction because of the interfacial compliance. The 
stiffness from the substrate alone (filled symbols) scales with 
radius, as expected from Hertz theory. Numerical uncertain- 
ties are comparable to the symbol size. 



B. Variations with load 

Figure El shows the load dependence of (a) normal dis- 
placement, (b) radius, (c) friction and (d) lateral stiffness 
for the same tips as Fig. (21 Each quantity is raised to 
a power that is chosen so that Hertz theory predicts the 
indicated straight line. A small finite-depth correction 
(~ 2%) is applied to the Hertz prediction for S (Eq. ||Sll). 

As also found for cylindrical contacts ^3 1 the normal 
displacement shows the smallest deviation from Hertz 
theory because it represents a mean response of many 
substrate atoms. Results for all bent crystals are nearly 
indistinguishable from the straight line. Results for the 
stepped surface are lower at small loads. Since the entire 
tip bottom contacts simultaneously, it takes a larger load 
to push the tip into the substrate. The amorphous results 
are shifted upwards by a fairly constant distance of about 
0.2 a. We define the origin of S as the tip height where 



the first substrate atom exerts a repulsive force on the 
tip. This is strongly dependent on the height of the low- 
est tip atom, while subsequent deformation is controlled 
by the mean tip surface. Agreement with Hertz is greatly 
improved by shifting the amorphous curve by this small 
height difference. Note that the zero of S is difficult to 
determine experimentally and is usually taken as a fit pa- 
rameter. If this is done, even results for the amorphous 
system yield values of R and E* that are within 10% of 
their true values. Thus measurements of 6, interpreted 
with continuum theory for spheres, can provide useful 
values of elastic properties at nanometer scales. 

As expected from the observed pressure distributions 
(Fig. 121), the contact radius is generally larger than the 
Hertz prediction. The shift is smallest for the dense tip 
because it approximates a continuous surface and the 
high density leads to a repulsive potential that rises more 
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than a hundred times more rapidly than for other tips. 
Results for bent crystal and amorphous tips are shifted 
upwards by fairly load-independent offsets of ^ 1 — 3 ti, 
leading to large fractional errors at low loads (up to 
100%). The stepped crystal shows qualitatively differ- 
ent behavior, with a rising in discrete steps as sequential 
terraces come into contact. Note that the size of the 
first terrace is not unique, but depends on the registry 
between the bounding sphere and crystalline lattice [4g . 
Larger deviations may be observed when the first step 
has very few atoms. Such tips may be more likely to be 
chosen for AFM studies because they tend to give sharper 
images. 



In order to predict the friction between surfaces, one 
must make assumptions about how F depends on area 
and load. The straight line in Fig. |2tc) corresponds to 
a friction force that is proportional to load. Static fric- 
tion values for bent and stepped commensurate surfaces 
are consistent with this line and a coefhcient of friction 
/i = F/N = 0.63. Analytic [23 and simple numerical |^ 
models show that this is a general feature of mated sur- 
faces where each tip atom has the same lateral position 
relative to substrate atoms. The friction on amorphous 
and incommensurate surfaces is always lower and scales 
more closely with the contact area, as indicated by broken 
line fits to oc N"^/^ and discussed further in Sec. IV BI 
|59| | . Many authors have made this assumption in fitting 
experimental data, but it is not obvious why it should 
hold. The friction per unit area between flat amorphous 
surfaces decreases as the square root of the area, but 
rises linearly with the normal pressure js^ ]. Wenning 
and Miiser have noted that these two factors may com- 
bine for spherical contacts to give a net friction that rises 
linearly with area [3ll |. However, their argument would 
predict that the frictional force in a cylinder-on-flat ge- 
ometry would not scale with area, and our previous simu- 
lations found that it did . Continuum theory predicts 
that the lateral stiffness k = 8G*a, and should follow the 
straight line in Fig. Ofd). Measured values of the total 
stiffness (open symbols) are always substantially lower. 
This is because continuum theory assumes that there is 
no displacement at the interface, only within the sub- 
strate. In reality, the frictional force is always associated 
with atomic scale displacements of interfacial atoms rel- 
ative to a local energy minimum [sol Is^ . The derivative 
of force with displacement corresponds to an interfacial 
stiffness ki that adds in series with the substrate contribu- 
tion (Eq. O) [Hi^. Our numerical results show that 
ki can reduce k by more than an order of magnitude, 
particularly for tips where F is small. We also evalu- 
ated the substrate stiffness fcsub by forcing all contacting 
substrate atoms to move together. These results (filled 
symbols) lie much closer to the Hertz prediction. Only 
the stepped tip shows large discrepancies, and these are 
correlated with the large deviation between the measured 
and predicted contact radii. 
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FIG. 4: (Color online) Local pressure vs. radius for bent 
commensurate tips with rcut = 1.5 a (filled squares) and 
Tcut = 2.2 cr (open squares) and with R — 100 a, ei/e = 0.5 
and N/{R^E*) = 0.0016. Dashed lines show Maugis-Dugdale 
fits with the measured work of adhesion and ho fit to the range 
of the interactions. A dotted line shows a fit for rcut = 2.2 a 
with /lo = 0.8 (T that improves agreement with the measured 
pulloff force (Fig. |HJ. Characteristic radii Ta, Vb and are 
defined by the locations where the radially averaged pressure 
first becomes zero, is most negative, and finally vanishes, re- 
spectively. 



V. ADHESIVE CONTACTS 
A. Pressure distribution 

Figure 0] compares the calculated pressure distribution 
in adhesive contacts with the Maugis-Dugdale prediction 
(lines). A bent commensurate tip was used to minimize 
deviations from continuum predictions for a sphere. Re- 
sults for two different rcut are presented to indicate diffi- 
culties in fitting longer-range interactions to M-D theory. 
The work of adhesion was calculated for unbent surfaces 
(Sec. IIII|I with Ei/e — 0.5, yielding w — 1.05 e/cr^ and 
1.65 e/cr^ for rcut = 1-5 a and 2.2 a, respectively. This 
leaves only one fitting parameter in M-D theory. For 
the dashed lines, the width of the attractive interaction 
ho = w/cTo was chosen to coincide with the range of 
the atomic potential. The dotted line shows a fit with 
ho = 0.8 a for r^ut = 2.2 cr, which gives better values for 
the pulloff force, but poorer radii (Sec. IVB|I . 

In M-D theory, it is common to identify two radii, a 
and c, with the inner and outer edges of the plateau in 
the pressure, respectively . For r < a the surfaces 

are in hard-sphere contact, and for a < r < c they are 
separated and feel the constant attraction. The contin- 
uously varying interactions between atoms in our sim- 
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ulations lead to qualitatively different behavior. There 
is no sharp transition where the surfaces begin to sepa- 
rate, and the attraction shows a smooth rise to a peak, 
followed by a decay. To facilitate comparison to contin- 
uum theories, we introduce the three characteristic radii 
indicated by arrows for each rcut- The innermost, Tq, cor- 
responds to the point where the interaction changes from 
repulsive to attractive and can be calculated in M-D the- 
ory. The outermost, rc, corresponds to c in M-D theory 
- the point where interactions end. The middle, r^, cor- 
responds to the peak attraction. It was also studied by 
Greenwood [S^ for contact of a featureless sphere that 
interacted with the flat substrate through Lennard-Jones 
interactions. He found n lay close to the JKR prediction 
for contact radius at large loads. All three radii converge 
in the limit of repulsive interactions or the JKR limit of 
an infinitely narrow interaction range. 

At small radii the atomistic results for p lie above M- 
D theory, and they drop below near r/a = 10. These 
deviations can be attributed to the increasing stiffness 
of the repulsion between tip and substrate atoms with 
increasing pressure. Just as in the non-adhesive case, the 
stiffer interactions in the center of the tip lead to bigger 
pressures for a given interpenetration. The change in 
pressure with separation produces less smearing at the 
edge of the repulsive region of the contact than for the 
non-adhesive case (Fig. El^b)), and the values of ra are 
very close to M-D theory. This can be understood from 
the fact that surfaces separate much more rapidly in JKR 
theory (oc (r — a)^^^) than in Hertz theory (ex (r — a)^^^) 
(Sec. II V|) . Thus the same height change in the repulsive 
region corresponds to a much smaller change in radius. 
Of course the finite range of the attractive tail of the 
potential leads to a broad region of adhesive forces out to 
Tc- The continuous variation of p in the attractive tail is 
only crudely approximated by M-D theory. The difficulty 
in determining the optimum choice for /iq increases with 
the range of interactions, as discussed further below. 

Figure |S1 shows the effect of tip geometry on pressure 
distribution. We found that the work of adhesion was 
very sensitive to tip geometry. To compensate for this ef- 
fect, we varied (Table I) to yield the same w for the tips 
in Fig. [S] Then all tips should have the same pressure 
distribution in continuum theory. The M-D predictions 
for p are not shown because even the bent commensurate 
tips produce significantly different results (Fig. 2J|. In- 
stead, we compare other tips to the bent commensurate 
tip. 

As for non- adhesive tips, local fiuctuations in pressure 
are small for commensurate tips (Fig. [S{a) and (d)) and 
comparable to the mean pressure for incommensurate or 
amorphous tips (Fig. Elb) and (c)). Note however, that 
the fluctuations become smaller in the adhesive regime 
(large r) . This is because the potential varies more slowly 
with height, so fluctuations in the separation between 
atoms have less effect on p. One consequence is that the 
outer edge of the contact is nearly the same for com- 
mensurate and incommensurate tips. The radii for the 



TABLE I: Relations between interaction strength e^, work of 
adhesion w, Nc and dimensionless pulloff force Nc/ ijiwR) for 
tips with different atomic-scale geometries. The last column 
gives the dimensionless pulloff force from M-D theory, A^^^~^, 
with the measured w and ho — 0.5 a. Values are accurate to 
the last significant digit. 



Tip geometry 












Commensurate 


0.213 


0.46 


256 


1.77 


1.74 




0.5 


1.05 


569 


1.72 


1.64 


Incommensurate 


0.535 


0.46 


258 


1.79 


1.74 




0.5 


0.45 


238 


1.68 


1.74 


Amorphous 


1.0 


0.46 


326 


2.26 


1.74 




0.5 


0.23 


136 


1.88 


1.83 


Stepped 


0.213 


0.46 


104 


0.72 


1.74 




0.5 


1.05 


168 


0.51 


1.64 



amorphous tip are significantly larger than bent tips, pre- 
sumably because of a much greater effective roughness. 
As for nonadhesive tips, the mean pressure on incom- 
mensurate tips is close to the commensurate results. Ad- 
hesion reduces the roughness-induced drop in pressure in 
the central region of the amorphous tip. For the stepped 
tip, the contact radius is dominated by the size of the 
terraces, but adhesive tails are visible outside the edge 
of each terrace. Only rc is easily defined for the stepped 
tips. This increases in a nearly stepwise manner, and its 
load dependence is not shown below. For the amorphous 
and incommensurate tips, radii are determined from the 
mean pressure at a given radius (open circles). Errorbars 
are less than 0.5 ct. 



B. Variations of radius and displacement with load 

Figure El^b) compares the measured contact radii for 
the tips of Fig. 13 to M-D theory as load is varied. The 
simulation results for (open symbols) decrease with 
decreasing load as (closed symbols) decreases to zero. 
All interactions in the contact are then adhesive. As 
Tc continues to drop, the area contributing to adhesion 
drops, and the load rises back toward zero. This regime is 
not considered in the original M-D theory. The extension 
to Tq = by Johnson and Greenwood (JG) is shown 
by a dashed line in the figure. Along this line the stress 
in the contact is constant, giving TV — —Uf^irr^. 

As for non-adhesive tips, the contacts tend to be larger 
than continuum theory predicts. However, the shift for 
bent tips is smaller than in the non-adhesive case, and the 
commensurate and incommensurate results are closer, as 
expected from Fig. [S] Larger deviations are observed for 
the amorphous tip, with radii typically 2 or 3 cr larger 
than predicted. The deviation becomes most pronounced 
at negative loads, where the amorphous tip remains in 
contact well below the predicted pulloff force. 

Figure EJa) compares the value of to the JKR pre- 
diction for contact radius. As found by Greenwood [52j |. 
the numerical results are close to JKR at large loads, but 
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FIG. 5: (Color online) Normal pressure vs. radius for four 
different tip geometries: (a) bent commensurate crystal, (b) 
bent incommensurate crystal, (c) amorphous, and (d) stepped 
crystal. In all cases the tip radius R — 100 a, the nor- 
malized load N / {B? E") = 0.0023 and the surface energy 
w — 0.46 e/(T^. Dots show the normal pressure on each sur- 
face atom. In (b) and (c), circles show the mean pressure in 
radial bins of width a, and lines show the bent commensurate 
results from (a). Horizontal dashed lines are at zero pressure. 



deviate at negative loads because M-D and JKR predict 
different puUofF forces. Since JKR assumes a singular ad- 
hesive stress at the radius of the contact, it seems natural 
that its predictions lie closest to the position of the peak 
tensile stress. 

Figure 13 shows the calculated radii for bent and amor- 
phous tips with the same interaction energy e^/e — 0.5. 
The small changes in tip geometry lead to a roughly four 
fold variation in both w and Nc (Table I). The largest w 
is obtained for commensurate tips because all atoms can 
optimize their binding coherently. Atoms on incommen- 
surate tips sample all positions relative to the substrate, 
and can not simultaneously optimize the binding energy. 
The larger height fluctuations on amorphous tips lead to 
even greater reductions in w. 

In the simulations, the pulloff force, iVc, corresponds 
to the most negative load where the surfaces remain in 
contact and the can be defined. Its normalized value, 
Nc/nwR, is equal to -3/2 in JKR theory, -2 in DMT the- 
ory, and lies between these values for M-D theory. Table 
I shows normalized results for various tips. As expected 
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FIG. 6: Variation of contact radii with external load for 
bent commensurate (squares), bent incommensurate (trian- 
gles) and amorphous (circles) tips with R — 100 a. In (a), 
the radius rj, where the force is most attractive is compared to 
JKR theory (solid line). In (b), values of Va (filled symbols) 
and Tc (open symbols) are compared to M-D theory (solid 
lines) with the JG extension for when Va = (dashed line). 
All continuum fits use the independently measured surface en- 
ergy w = 0.46 e/cr^ and an interaction range ho = 0.5 a that 
is consistent with the potential range. Numerical uncertainty 
in radii is comparable to the symbol size. 



from the good agreement in Figs. and [7| results for 
bent tips lie between JKR and DMT values and can be 
fit with M-D theory. The values for stepped and amor- 
phous tips lie outside the bounds of M-D theory. This is 
an important result because pulloff forces are frequently 
used to measure the work of adhesion. Based on con- 
tinuum theory, one would expect that the uncertainty in 
such measurements is less than the difference between 
JKR and DMT predictions. Our results show factor of 
two deviations for stepped tips, which may be common in 
scanning probe experiments. Other simulations showed 
that the stepped tip values are strongly dependent on the 
size of the first terrace, as well as any tilt of the terraces 
or incommensurability. 

It might seem surprising that the stepped tip has a 
smaller pulloff force than the bent tip, because the en- 
tire first terrace can adhere without any elastic energy 
penalty. However this effect is overcome by the greater 
contact size for bent tips: The radius of the first terrace, 
Tt ~ 6 (7, is smaller than the values of and Vc at pulloff 
for bent tips. As the adhesion is decreased, the predicted 
contact size at pulloff will drop below rt and the stepped 
tip may then have a larger pulloff force than predicted. 
This limit can also be reached by increasing the width 
of the first terrace. For a tip that lies entirely within a 
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FIG. 7; Contact radii as a function of normal load for (a) 
bent commensurate, (b) bent incommensurate, and (c) amor- 
phous tips with Ei = 0.5e. While all have the same interaction 
strength, the work of adhesion varies by more than a factor 
of four. Broken lines show M-D predictions for (triangles) 
and Tc (circles) with the indicated values of w and ho — 0.5 a. 
The value of rt (squares) is closest to the JKR result (solid 
line). Deviations from the continuum predictions show the 
same trends as in Fig. |S| Numerical uncertainty in radii is 
comparable to the symbol size. 



sphere of radius R, r'f < — {R — ~ 2dR where d 
is the layer spacing in the crystal. For our geometry this 
corresponds to rj < 12 cr, which is about twice the value 
for our tip. As noted above, terraces with even smaller 
radii may be preferentially selected for imaging studies 
and will lead to lower |7Vc|. 

As Tcut increases, it becomes harder for the M-D model 
to simultaneously fit both the radii and the puUoff force. 
Figure IHl shows simulation data for a bent commensurate 
tip with Tcut — 2.2 cr. Using the value of ho — 1.2 cr (Fig. 
^ reproduces all radii fairly well at large loads, but gives 
a substantially different pulloff force, —906 e/a instead 
of —872 e/(T. Decreasing ho to 0.8 a fits the pulloff force, 
and improves the fit to ra at low TV. However, the pre- 
dicted values for Tq at large N are slightly too high and 
the values for Tc are shifted far below (~ 2 — 3 cr) simula- 
tion data. This failure is not surprising given the crude 
approximation for adhesive interactions in M-D theory. 
As might be expected, the best values for the pulloff force 
are obtained by fitting the region near the peak in the 
force, rather than the weak tail (Fig. ^J. 

It should be noted that for bent crystalline and amor- 
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FIG. 8: Fits of M-D theory (broken lines) to Tc (squares) and 
Va (circles), and of JKR theory (solid line) to Vb (triangles) for 
a bent commensurate tip with rcut = 2.2 a. For ho — 1.2 a 
(dashed line) M-D theory fits both and Va at large loads, 
but the magnitude of the pulloff force is too large. Decreasing 
ho to 0.8 a (dotted lines) gives excellent agreement with the 
pulloff force, but not the radii (Fig. |1J. Numerical uncertainty 
in radii is comparable to the symbol size. 



phous tips, all of our results for r can be fit accurately 
to M-D theory if E*, w, and R are taken as adjustable 
parameters. The typical magnitude of adjustments is 10 
to 30%, which is comparable to typical experimental un- 
certainties for these quantities in nanoscale experiments. 
Indeed one of the common goals of these experiments 
is to determine any scale dependence in continuum pa- 
rameters. For this reason it would be difhcult to test 
continuum theory in scanning probe experiments even if 
r could be measured directly. 

Experiments can more easily access the variation of 
normal displacement with load. This requires subtract- 
ing the height change due to the normal compliance of 
the machine controlling the tip, which is difficult for 
standard AFM cantilevers, but possible in stiffer devices 
fl^ Il3| . The absolute zero of 6 is not fixed, but must 
be fitted to theory. Figure shows two measures of the 
normal displacement in our simulations. One, S^p, corre- 
sponds to the experimentally accessible tip displacement. 
We associate the zero of (Stip with the point where the first 
substrate atom exerts a repulsive force. The second, (5sur, 
is the actual depression of the lowest substrate atom on 
the surface relative to the undeformed surface. The two 
differ because of the interfacial compliance normal to the 
surface, which is assumed to vanish in continuum theory. 

The simulation results for S are more sensitive to the fi- 
nite sample depth D and lateral periodicity L than other 
quantities. To account for this, the predictions of M-D 
and JKR theory are shifted by the leading analytic cor- 
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FIG. 9; Normal displacement as a function of external load 
measured from (a) the depression of the lowest substrate atom 
5sur and (b) the displacement of the tip relative to the height 
where the first substrate atom exerts a repulsive force Stip- 
The tips are the same as in Fig. |5] The JKR prediction is in- 
dicated by dotted lines, the M-D prediction by solid lines, and 
the JG extension for Va = with dashed lines. All are cor- 
rected for the finite dimensions of the substrate as described 
in the text, and uncertainties are comparable to the symbol 
size. The difference AS = Stip — 5sur reflects the interfacial 
compliance, and is shown for the bent incommensurate tip in 
the inset. 



rections in a/D |37ll3ll39l|: 

6 = Sh(1 - b^) + (5adhcsion(l - d^) , (8) 

where 6h is the Hertz prediction (Eq. Q), and d are 
fit parameters, and (^adhesion — —^2wTra/E* for JKR 
theory and —{2ao/E*)\/c'^ — for M-D theory. We ob- 
tained b — 0.8 from simulations with dense, non-adhesive 
tips (Fig. I2J) and d = 0.3 from simulations with dense, 
adhesive tips. Results for dense tips are then indistin- 
guishable from the fit lines in Figs. Of a) and El Values 
of b and d from numerical studies of continuum equa- 
tions are of the same order [sSEEESj but to our knowl- 
edge these calculations have not considered our geometry 
where L « Z?. 

With the finite-size corrections, continuum results for 6 
lie very close to the simulation data for bent tips. Agree- 
ment is best for Sgui- because neither it nor continuum 
theory include the interfacial compliance. The choice of 
zero for 5tip is not precisely consistent with M-D theory. 
The first repulsive interaction would correspond to the 
first non-zero , which occurs at a slightly negative S in 
M-D theory. For the parameters used here, the shift in S 
is about 0.24 a, which is about twice the discrepancy be- 
tween Stip and M-D theory at large loads. This implies 
that the interfacial compliance produces a comparable 



correction. 

The effects of interfacial compliance are biggest for the 
most negative values of i5tip. Here ^tip decreases monoton- 
ically as N increases to zero. In contrast, Ssnr is flat and 
then increases back towards zero. In this regime, ra — 
and the only interaction comes from the attractive tail 
of the potential. The net adhesive force gradually de- 
creases {N increases) as the magnitude of the separation 
between tip and substrate, AS = Sup — (5sur, increases 
(inset). Most of this change occurs in i5tip, while the 
displacement of the surface relaxes back to zero as the 
attractive force on it goes to zero. The JG extension to 
M-D theory is expressed in terms of Sup and provides a 
good description of its change with N (dashed line) , even 
with the assumption of a constant attractive (Tq. 

For amorphous tips, both values of S are substantially 
above M-D theory. The shifts are bigger than in the non- 
adhesive case, about 0.35 a vs. 0.25 a. This is correlated 
with the larger than predicted puUoff force. Based on the 
increase in \Nc\, the effective work of adhesion appears 
to be larger by about 30% than that measured for a flat 
surface. The stronger adhesive contribution to the total 
normal force leads to a larger value of S. One may un- 
derstand these changes in terms of the effect of surface 
roughness. Atomic-scale roughness on a nominally flat 
surface prevents many atoms from optimizing their bind- 
ing energy. As S and the contact area shrink, the long 
wavelength height fluctuations become irrelevant and no 
longer prevent the few atoms remaining in the contact 
from adhering. Thus while the large load values of S can 
be flt to continuum predictions with the measured w and 
a simple shift in the origin of S, the small N values cor- 
respond to a larger work of adhesion. The magnitude of 
the increase (~ 30%) is modest given that the incommen- 
surate tip has about twice as large a w as the amorphous 
tip for the same interaction energy e^. 

The data for stepped tips are qualitatively different 
than the others. As for the non-adhesive case, 5 is lower 
than the continuum prediction at large loads, because 
the flat tip is harder to push into the substrate. The de- 
viation increases rapidly at negative loads, with a sharp 
drop near = where the contact shrinks to the flrst 
terrace. As noted above, the radius of the first terrace is 
smaller than the radius at puUoff predicted by continuum 
theory, and the puUoff force is less than half the predicted 
value. 



C. Friction and lateral stiffness 

Scanning probe microscopes can most easily de- 
tect friction forces and the lateral compliance of the 
tip/substrate system 0, |^. Figure El shows F as a 
function of load for five different adhesive tips. As for 
non-adhesive tips, tip geometry has a much larger effect 
on F than other quantities, and values for bent commen- 
surate and incommensurate tips differ by two orders of 
magnitude. 
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FIG. 10: (Color online) Static friction F vs. load N for 
the indicated tip geometries with same systems as in Fig. |S] 
Attempts to fit the data by assuming F is proportional to the 
area predicted by JKR theory and the interpolation scheme 
of Ref. I29II are shown by dotted and solid lines, respectively. 
Dashed lines are separate linear fits for the stepped tip for 
the cases where one or two terraces are in repulsive contact. 
Numerical uncertainties are comparable to the symbol size. 



Since the friction was measured at constant load (Sec. 
only values in the stable regime, dSnp/dN > could 
be obtained. Even in this regime, we found that the tip 
tended to detach at TV > A'^c. This was particularly pro- 
nounced for commensurate tips. Indeed the bent and 
stepped commensurate tips detached after the first peak 
in the friction force for all the negative loads shown in 
Fig. Iiur a'). At loads closer to the puUoff force, detach- 
ment occurred at even smaller lateral displacements. 

As noted above, bent commensurate tips have the 
strongest adhesion energy because all atoms can simul- 
taneously optimize their binding. For the same reason, 
the adhesion energy changes rapidly as atoms are dis- 
placed laterally away from the optimum position, al- 
lowing puUoff above the expected Nc- The extent of 
the charige depends on the sliding direction and registry 
Es"; so] . We consider shding in the (100) direction 
(Sec. imp , where atoms move from points centered be- 
tween substrate atoms towards points directly over sub- 
strate atoms. This greatly reduces the binding energy, 
leading to detachment at less than half the puUoff force. 
Note that changes in binding energy with lateral displace- 
ment lead directly to a lateral friction force [s^l and the 
bent and stepped commensurate tips also have the high- 
est friction. We suspect that tips with higher friction 
may generally have a tendency to detach farther above 



TABLE II: Frictional stresses Tf and apparent moduli G" eval- 
uated from the derivatives of fits in Figs. Illl and lT^ respec- 
tively, as a function of tip geometry and work of adhesion w. 
Values of G' are less than the effective shear modulus of the 
substrate, G* = 18.3 e/a^, and frictional stresses are much 
smaller than expected from continuum expressions (Fig. 1131 . 
Statistical errorbars are indicated in parentheses. 



Tip geometry 


w{ea ^) 




G'iea--) 


Commensurate 





1.35(7) 


13.7(4) 




0.46 


1.82(4) 


12.8(3) 


Incommensurate 





0.0044(7) 


0.70(5) 




0.46 


0.0151(8) 


1.0(3) 


Amorphous 





0.056(4) 


1.53(5) 




0.46 


0.24(2) 


4.3(6) 



Nc during sliding than other tips. 

At the macro scale, friction is usually assumed to be di- 
rectly proportional to load. All the tips have substantial 
friction at zero load, due to adhesion. The friction force 
also varies non-linearly with N, showing discrete jumps 
for the stepped tip, and curvature for the other tips, par- 
ticularly near Nc- The curvature at small N is reminis- 
cent of the dependence of radius on load. Several authors 
have fit AFM friction data assuming that F scales with 
contact area, and using a continuum model to determine 
the area as a function of load HlllElllililliini- The 
dotted lines in Fig. 1101 show attempts to fit F using the 
area from JKR theory with w adjusted to fit the puUoff 
force and the proportionality constant chosen to fit the 
friction at high loads. None of the data is well fit by 
this approach. The solid lines show fits to an expression 
suggested by Schwarz 29] that allows one to simply in- 
terpolate between DMT and JKR limits as a parameter 
ti is increased from to 1. Reasonable fits are possible 
for non-stepped tips when this extra degree of freedom is 
allowed. Values of the friction per unit area are within 
25% of those determined below from the true area (Ta- 
ble II), but the values of w and ti do not correspond to 
any real microscopic parameters. For example, for the 
three non-stepped tips where a direct measurement gives 
w = 0.46 e/a^, the fits give w = 0.41 e/a^, 0.55 e/a^ 
and 0.94 e/tr^ for commensurate, amorphous and incom- 
mensurate tips, respectively. The value of ti varies from 
0.1 to 0.5 to 1.5, respectively. Note that the value of 
ti = 1.5 for incommensurate tips is outside the physi- 
cal range. In this case, and some experiments the 
friction rises more slowly than the JKR predictioii, while 
M-D theory and simpler interpolation schemes |2^ Isof 
always give a steeper rise. Such data seem inconsistent 
with F scaling with area. 

Our simulations allow us to test the relationship be- 
tween F and area without any fit parameters. However, 
it is not obvious which radius should be used to deter- 
mine area. Figure 1111 shows friction plotted against , 

and r^. The stepped tip is not shown, since only Vc is 
easily defined and it increases in one discontinuous jump 
over the range studied. For all other tips the friction 
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FIG. 11: (Color online) Static friction F plotted against 
radius squared for the indicated tip geometries with same pa- 
rameters as in Fig. |S| Values of (filled circles), rl (open 
squares), and (filled triangles) are shown for adhesive tips. 
Open circles show F vs. for non-adhesive tips with data in 
(b) and (c) multiplied by a factor of two for clarity. Dashed 
lines show unconstrained linear fits to each data set. Numer- 
ical uncertainties are comparable to the symbol size. 



is remarkably linear when plotted against any choice of 
radius squared. In contrast, plots of N vs. show sig- 
nificant curvature. 

For F to be proportional to area, the curves should 
also pass through the origin. This condition is most 
closely met by r^, except for the incommensurate case. 
The idea that friction should be proportional to the area 
where atoms are pushed into repulsion seems natural. 
The other radii include attractive regions where the sur- 
faces may separate far enough that the variation of force 
with lateral displacement is greatly reduced or may even 
change phase. However, in some cases the extrapolated 
friction remains finite as goes to zero and in others it 
appears to vanish at finite r^. 

Also shown in Fig. ^2 are results for non-adhesive 
tips (open circles) . Values for incommensurate and amor- 
phous tips were doubled to make the linearity of the data 
more apparent. Only the bent commensurate data shows 
significant curvature, primarily at small N . As noted in 
Sec. IIVI friction is proportional to load for this tip, and 
the curvature is consistent with this and ~ N"^^^. No 
curvature is visible when adhesion is added, but the fact 
that the linear fit extrapolates to F = at ra/a ~ 4 sug- 
gests that the linearity might break down if data could 



be obtained at lower loads. 

Results for adhesive and non-adhesive amorphous tips 
can be fit by lines through the origin (not shown), 
within numerical uncertainty. The same applies for non- 
adhesive bent incommensurate tips, but adding adhesion 
shifts the intercept to a positive force at Va — 0. It would 
be interesting to determine whether this extrapolation is 
valid. Friction can be observed with purely adhesive in- 
teractions, since the only requirement is that the magni- 
tude of the energy varies with lateral displacement (50j. 
However, it may also be that the friction on incommen- 
surate tips curves rapidly to zero as ra — > 0. Indeed one 
might expect that for very small contacts the tip should 
behave more like a commensurate tip, leading to a more 
rapid change of F with area or load. The only way to ac- 
cess smaller ra is to control the tip height instead of the 
normal load. This is known to affect the measured fric- 
tion force "4^ , and most experiments are not stiff enough 
to access this regime. 

FigurelT^shows the lateral stiffness as a function of the 
three characteristic radii. Except for the stepped tip (not 
shown), k rises linearly with each of the radii. As for non- 
adhesive tips, the slope is much smaller than the value 
8G* predicted by continuum theory (solid lines) because 
of the interfacial compliance (Eq. {T))). The intercept is 
also generally different from the origin, although it comes 
closest to the origin for ra- As for friction, it seems that 
the repulsive regions produce the dominant contribution 
to the stiffness. 

Results for non-adhesive tips are also included in Fig. 
[T^ They also are linear over the whole range, and the fits 
reach A: = at a finite radius Ofc. This lack of any stiff- 
ness between contacting surfaces seems surprising. Note 
however that linear fits to Fig. |SJb) also would suggest 
that the radius approached a non-zero value, ao, in the 
limit of zero load. Moreover, the values of oq and ak fol- 
low the same trends with tip geometry and have similar 
sizes. 

The finite values of oq and Uk can be understood from 
the finite range of the repulsive part of the interaction. 
As long as atoms are separated by less than Tcut there 
is a finite interaction and the atoms are considered in- 
side a. However, the force falls rapidly with separation, 
and atoms near this outer limit contribute little to the 
friction and stiffness. If 6h is the distance the separa- 
tion must decrease to get a significant interaction, then 
ao = {2R6hy/^ at the point where the first significant 
force is felt. Taking the estimate of Sh = 0.04 a from 
Sec. IIVI then aq ~ 3 cr, which is comparable to the ob- 
served values of aq and a^. The shift is smaller for the 
adhesive case because there are still strong interactions 
when Tq goes to zero. The larger shifts for the amor- 
phous tip may reflect roughness, since the first point to 
contact need not be at the origin. We conclude that the 
linear fits in Fig. ^Igo to zero at finite radius, because 
ra overestimates the size of the region that makes signifi- 
cant contributions to forces, particularly for non-adhesive 
tips. Note that the plots for friction with non-adhesive 
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FIG. 12: (Color online) Lateral stiffness k plotted against 
radius for the indicated tip geometries with same parameters 
as in Fig. |^ Values of Va (filled circles), (open squares), 
and Tc (filled triangles) are shown for adhesive tips. Open 
circles show F vs. a for non-adhesive tips. Broken lines show 
unconstrained linear fits to each data set, and the solid line 
indicates the slope predicted by continuum theory. Numerical 
uncertainties are comparable to the symbol size, except for the 
adhesive data in (b) where they may be up to 50%. 



tips (Fig. Ill|) are also consistent with an offset, but that 
the offset appears much smaller when plotted as radius 
squared. 

The slope of the curves in Fig. ^Jcan be used to define 
a differential friction force per unit area or yield stress 
Tf = dF/dirr'^. It is interesting to compare the magni- 
tude of these values (Table II) to the bulk yield stress of 
the substrate Ty. Assuming Lennard- Jones interactions, 
the ideal yield stress of an fee crystal in the same shear- 
ing direction is 4 to 10 ea^, depending on whether the 
normal load or volume is held fixed. The commensurate 
tip is closest to a continuation of the sample, and the 
force on all atoms adds coherently. As a result Tf is of 
the same order as Ty, even for the non-adhesive tip. Val- 
ues for adhesive amorphous and incommensurate tips are 
about one and two orders of magnitude smaller than Ty , 
respectively. This reflects the fact that the tip atoms can 
not optimize their registry with the substrate. Removing 
adhesive interactions reduces Ty by an additional factor 
of about four in both cases. 

In continuum theory, k — 8G*r, and the slope of fits 
in Fig. El could then be used to determine the effective 
shear modulus G*. However, as noted above, the in- 
terfacial compliance leads to much lower stiffnesses. To 
illustrate the magnitude of the change we quote values 
of G' = {l/8)dk/dr in Table II. AU values are below the 



true shear modulus G* = 18.3 e/cr^ obtained from the 
substrate compliance alone (Fig. O. As always, results 
for bent commensurate tips come closest to continuum 
theory with G"/G* ~ 0.7. Values for adhesive amorphous 
and incommensurate tips are depressed by factors of 3 
and 20 respectively, and removing adhesion suppresses 
the value for amorphous tips by another factor of four. 

Carpick et al. noted that if friction scales with area 
and k with radius, then the ratio F/k'^ should be con- 
stant 0, IE 0j 113 ■ Defining the frictional force per unit 
area as t/ and using the expression for k from continuum 
theory, one finds 6AG*^F/-Kk^ = Tf - principal, this al- 
lows continuum predictions to be checked and r/ to be 
determined without direct measurement of contact size. 
Figure [T^ shows: 



^cff 



64(G*)^F/7rfc' 



(9) 



as a function of for different tip geometries and in- 
teractions. Except for the non-adhesive stepped tip, the 
value of F/k'^ is fairly constant at large loads, within 
our numerical accuracy. Some of the curves rise at small 
loads because the radius at which F reaches zero in Fig. 
^2 tends to be smaller than that where k reaches zero in 
Fig. 1121 These small radii are where continuum theory 
would be expected to be least accurate. Note that the 
deviations are larger for the non-adhesive tips, perhaps 
because the data extends to smaller radii. 

The data for stepped tips are of particular interest be- 
cause the contact radius jumps in one discrete step from 
the radius of the first terrace to the radius of the second. 
The friction and stiffness also show discontinuous jumps. 
Nonetheless, the ratio F/k^ varies rather smoothly and 
even has numerical values close to those for other tips. 
The most noticeable difference is that the data for the 
nonadhesive stepped tip rises linearly with load, while 
all other tips tend to a constant at high load. These re- 
sults clearly demonstrate that success at fitting derived 
quantities like F and k need not imply that the true con- 
tact area is following continuum theory. 

The curves for Ty* in Fig. 1 131 are all much higher than 
values of the frictional stress Tf obtained directly from 
the friction and area (Table II). Even the trends with 
tip structure are different. The directly measured fric- 
tional stress decreases from bent commensurate to amor- 
phous to bent incommensurate, while Tf^ is largest for 
the amorphous and smallest for the bent commensurate 
tip. These deviations from the continuum relation are 
directly related to the interfacial compliance ki . The con- 
tinuum expression for the lateral stiffness neglects ki and 
gives too small a radius at each load. This in turn over- 
estimates the frictional stress by up to two orders of mag- 
nitude. Similar effects are likely to occur in experimental 
data. 

Experimental plots of F/fc^ have been obtained for 
silicon- nitride tips on mica and sodium chloride [E 0, 113 
and on carbon fibers j^]- Data for carbon fibers and 
mica in air showed a rapid rise with decreasing N at 
low loads M, I6Q| ■ For mica the increase is almost an or- 
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FIG. 13: (Color online) Ratio of friction to stiffness squared 
as a function of load for the indicated tip geometries (a) with- 
out and (b) with adhesion (w = 0.46 e/o"^). Lines are guides 
to the eye. Numerical uncertainties are comparable to the 
symbol size, except for the bent incommensurate data in (b), 
where they may be as large as 50%. 



der of magnitude, which is comparable to our results for 
non-adhesive bent incommensurate tips. This correspon- 
dence may seem surprising given that the experiments 
measured adhesion in air. However, the adhesive force 
was mainly from long-range capillary forces that operate 
outside the area of contact. Following DMT theory, they 
can be treated as a simple additive load that does not 
affect the contact geometry. In contrast, data for mica 
in vacuum is well fit by JKR theory, implying a strong 
adhesion within the contact 0, Q • The measured value 
of F/k'^ is nearly constant for this system, just as in our 
results for most adhesive tips. Results for carbon fibers 
in vacuum show a linear rise like that seen for nonad- 
hesive stepped tips. 

From a continuum analysis of the carbon fiber data, 
the frictional stress was estimated to be t/ ~ 300 MPa 
assuming a bulk shear stress of G* = 9.5 GPa @. Note 
that Fig. US would suggest t//G* - 0.1 to 0.3, while 
the true values (Table II) are as low as 0.0002. The data 
on carbon fibers could be fit with the bulk shear mod- 
ulus, but data on mica and NaCl ^, ,6QJ indicated that 
G* was 3 to 6 times smaller than bulk values. Our re- 
sults show that the interfacial compliance can easily lead 
to reductions of this magnitude and a corresponding in- 
crease in Tf, and that care must be taken in interpreting 
experiments with continuum models. 



VI. DISCUSSION AND CONCLUSION 

The results described above show that many different 
effects can lead to deviations between atomistic behav- 
ior and continuum theory and quantify how they depend 
on tip geometry for simple interaction potentials (Fig. 

In general, the smallest deviations are observed for 
the idealized model of a dense tip whose atoms form a 
nearly continuous sphere, although this tip has nearly 
zero friction and lateral stiffness. Deviations increase as 
the geometry is varied from a bent commensurate to a 
bent incommensurate to an amorphous tip, and stepped 
tips exhibit qualitatively different behavior. Tip geom- 
etry has the smallest effect on the normal displacement 
and normal stiffness (Fig. |31and|21) because they reflect 
an average response of the entire contact. Friction and 
lateral stiffness are most affected (Fig. O and IIOII , be- 
cause they depend on the detailed lateral interlocking of 
atoms at the interface. 

One difference between simulations and continuum 
theory is that the interface has a finite normal compli- 
ance. Any realistic interactions lead to a gradual in- 
crease in repulsion with separation rather than an ide- 
alized hard-wall interaction. In our simulations the effec- 
tive range over which interactions increase is only about 
4% of the atomic spacing, yet it impacts results in several 
ways. For bent commensurate tips it leads to an increase 
in pressure in the center of the contact (Figs. |21and0J. 
The pressure at the edge of nonadhesive contacts drops 
linearly over about 2 a, while continuum theory predicts 
a diverging slope. The width of this smearing grows as 
the square root of the tip radius and leads to qualitative 
chaiiges in the probability distribution of local pressures 
EE EE 113 ■ These effects could be studied in con- 
tinuum theories with soft-wall interactions. The normal 
interfacial compliance also leads to the offset in linear fits 
of F vs. and k vs. (Figs. ITTI and [T^ . Fits to the 
friction and local stiffness extrapolate to zero at finite 
values of Va because atoms at the outer edge of the re- 
pulsive range contribute to Tq but interact too weakly to 
contribute substantially to F and k. This effect is largest 
for non-adhesive tips. 

Approximating a spherical surface by discrete atoms 
necessarily introduces some surface roughness. Even 
bent crystalline tips have atomic scale corrugations, re- 
flecting the variation in interaction as tip atoms move 
from sites nestled between substrate atoms to sites di- 
rectly above. Amorphous and stepped tips have longer 
wavelength roughness associated with their random or 
layered structures respectively. This longer wavelength 
roughness has a greater effect on the contacts. For non- 
adhesive interactions, incommensurate and amorphous 
tips have a lower central pressure and wider contact ra- 
dius than predicted for ideal spheres. These changes are 
qualitatively consistent with continuum calculations for 
spheres with random surface roughness 16] . However the 
effective magnitude of the rms roughness A is smaller 
than expected from the atomic positions. The correlated 
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deviations from a sphere on stepped tips, lead to qual- 
itative changes in the pressure distribution on the sur- 
face (Fig. [3 and El- However, these changes are also 
qualitatively consistent with what continuum mechanics 
would predict for the true tip geometry, which is closer 
to a flat punch than a sphere. We conclude that the 
usual approximation of characterizing tips by a single 
spherical radius is likely to lead to substantial errors in 
calculated properties. Including the true tip geometry 
in continuum calculations would improve their ability to 
describe nanometer scale behavior. Unfortunately this is 
rarely done, and the atomic-scale tip geometry is rarely 
measured. Recent studies of larger tips and larger scale 
roughness are an interesting step in this direction . 

Roughness also has a strong influence on the work 
of adhesion w (Table I). Values of w were determined 
independently from interactions between nominally flat 
surfaces. For a given interaction strength, commensu- 
rate surfaces have the highest w, because each atom can 
optimize its binding simultaneously. The mismatch of 
lattice constants in incommensurate geometries lowers 
w by a factor of two, and an additional factor of two 
drop is caused by the small (A ^ 0.3 a) height fluctua- 
tions on amorphous surfaces. In continuum theory, these 
changes in w should produce nearly proportional changes 
in puUoff force Nc, and tips with the same w and ho 
should have the same Nc- Measured values of Nc differ 
from these predictions by up to a factor of two. It is par- 
ticularly significant that the dimcnsionlcss puUoff force 
for amorphous and stepped tips lies outside the limits 
provided by JKR and DMT theory. Experimentalists of- 
ten assume that these bounds place tight limits on errors 
in inferred values of w. 

In the case of amorphous tips the magnitude of Nc is 
30% higher than expected. The higher than expected 
adhesion in small contacts may reflect a decrease in ef- 
fective roughness because long- wavelength height fluctu- 
ations are suppressed. Stepped tips show even larger 
deviations from continuum theory that are strongly de- 
pendent on the size of the first terraces . Tips selected 
for imaging are likely to have the smallest terraces and 
the largest deviations from continuum theory. 

Adding adhesion introduces a substantial width to the 
edge of the contact, ranging from the point where inter- 
actions first become attractive Va to the outer limits of 
attractive interactions Vc (Fig. As the range of inter- 
actions increases, it becomes increasingly difficult to fit 
both these characteristic radii and the puUoff force with 
the simple M-D theory (Fig. |SJ). For short-range inter- 
actions, good fits are obtained with the measured w for 
bent tips. Data for amorphous tips can only be fit by 
increasing w, due to the reduction in effective roughness 
mentioned above (Fig. 01. For stepped tips the contact 
radius increases in discrete jumps as successive terraces 
contact the surface. 

The normal interfacial compliance leads to significant 
ambiguity in the definition of the normal displacement 
as a function of load (Fig. |5Jl. Continuum theory nor- 



mally includes only the substrate compliance, while ex- 
perimental measures of the total tip displacement Sup 
include the interfacial compliance. The substrate compli- 
ance was isolated by following the displacement of sub- 
strate atoms, (5sur and found to agree well with theory for 
bent tips in the repulsive regime. Johnson and Green- 
wood's extension of M-D theory ^2§\ includes the interfa- 
cial compliance in the attractive tail of the potential. It 
provides a good description of ^tip in the regime where 
ra = 0. Here 5 tip — Ssur increases to the interaction range 
ho. Results for amorphous tips show the greater adhe- 
sion noted above. Stepped tips follow continuum theory 
at large loads but are qualitatively different at negative 
loads. 

The most profound effects of tip geometry are seen 
in the lateral stiffness k and friction F, which vary by 
one and two orders of magnitude respectively. Contin- 
uum theories for k do not include the lateral interfacial 
compliance fcj. This adds in series with the substrate 
compliance fcgub (Eq. 10). Except for commensurate 
tips, fcj << fcsub and the interface dominates the total 
stiffness Experiments have also seen a substantial 

reduction in the expected lateral stiffness from this effect 

The friction on non-adhesive commensurate tips (bent 
or stepped) increases linearly with load, as frequently ob- 
served for macroscopic objects. In all other cases, is a 
nonlinear function of load. Our ability to directly mea- 
sure contact radii allowed us to show that F scales lin- 
early with contact area for incommensurate, amorphous 
and adhesive bent commensurate tips. These tips also 
show a linear scaling of k with radius. While these scal- 
ings held for any choice of radius, the linear fits are offset 
from the origin. It appears that the effective area con- 
tributing to friction and stiffness is often a little smaller 
than the area of repulsive interactions corresponding to 
Ta- As noted above, the offset from ra appears to corre- 
spond to the finite range over which repulsive forces rise 
at the interface. 

Experimental data |i,lllllllEll|l3,lill for friction 
and stiffness have been fit to continuum theory with the 
assumptions that F o(. r'^ and k cx r, but without the 
offsets seen in Figs. IlllandlT^ We showed that our data 
for bent and amorphous tips could be fit in this way (Fig. 
llOfl . but that the fit parameters did not correspond to 
directly measured values. This suggests that care should 
be taken in interpreting data in this manner. 

We also examined the ratio F/fc^. In continuum the- 
ory, this is related to the friction per unit area t'^^ 
through Eq. ©. Our results for F/k'^ (Fig. EJl show the 
range of behaviors observed in experiments, with a rela- 
tively constant value for adhesive cases, a rapid increase 
at low loads in some nonadhesive cases, and a linear rise 
for non-adhesive stepped tips. The directly measured 
values of r/ (Table II) are smaller than tJ^ by up to 
two orders of magnitude, and have qualitatively different 
trends with tip geometry. The difference is related to a 
reduction in the stiffness k due to interfacial compliance. 
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This reduces the inferred value of bulk shear modulus G' 
and increases the calculated contact area at any given 
area. We expect that experimental results for F/k"^ will 
produce similar overestimates of the true interfacial shear 
force. 

It remains unclear why F and k should follow the ob- 
served dependence on r^. Analytic arguments for clean, 
flat surfaces indicate that F is very sensitive to struc- 
ture, with the forces on commensurate, incommensurate 
and disordered surfaces scaling as different powers of area 
[3^ Isoj l . Only when glassy layers are introduced between 
the surfaces, does the friction scale in a universal man- 
ner [33,|63,ES0I- Wenning and Miiser jsj have argued 
that the friction on clean, amorphous tips rises linearly 
with area because of a cancellation of two factors, but 
have not considered k. Naively, one might expect that 
the length over which the force rises is a constant frac- 
tion of the lattice spacing and that k is proportional to 
F. However, the friction traces change with load and do 
not always drop to zero between successive peaks. We 
hope that our results will motivate further analytic stud- 
ies of this problem, and simulations with glassy films and 



more realistic potentials. 

While we have only considered single asperity contacts 
in this paper, it is likely that the results are relevant 
more broadly. Many experimental surfaces have random 
roughness on all scales that can be described by self-affine 
fractal scaling. Continuum models of contact between 
such surfaces show that the radius of most of contacts 
is comparable to the lower length scale cutoff in fractal 
scaling [5^ [s^l . This is typically less than a micrometer, 
suggesting that typical contacts have nanometer scale di- 
mensions where the effects considered here will be rele- 
vant. 
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